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PREFACE 

This  document  is  the  result  of  2  years  of  independent  research  and  study  in 
the  area  of  minimum  distance  theory  and  its  application  to  geodesic  measure¬ 
ment  and  surface  navigation.  The  purpose  is  to  develop  a  new  approach  to 
solving  the  problem  of  getting  from  point  a  to  point  b  on  well  behaved  sur¬ 
faces.  The  study  is  confined  to  the  use  of  curvilinear  coordinates  (u^,  u^, 

Ug )  and  to  surfaces  which  can  be  represented  by  a  function,  k,  such  that: 

Ug  =  k(u^,  u^) ,  (ie,  surfaces  where  one  coordinate  can  be  written  as  a 
function  of  the  remaining  two). 

In  rederiving  the  relations  which  determine  the  geodesics  or  curves  of  minimum 
distance,  we  shall  see  how  using  parametric  separation  of  variables,  along 
with  the  above  coordinate  variables  and  surface  restrictions,  greatly  simpli¬ 
fies  the  problem,  thus  giving  manageable  solutions. 

The  paper  addresses  the  development  of  solutions  not  just  from  the  point  of 
theoretical  interest,  but  in  terms  of  product  useability,  for  eventual  soft¬ 
ware  implementation.  By  product  useability,  I  mean  we  need  to  know  more  than 
just  the  relations  between  U2  and  Ug  for  an  arbitrary  minima  curve,  we  need 

explicit  relations  for  distance  traveled  and  direction  (measured  normally  on 
the  surface)  of  the  minima  curve.  Thus,  we  consider  a  closed  set  of  equations 
the  solution  for  a  given  surface  (and  initial  and/or  boundary  conditions)  to 
be  of  the  form: 

U3  =  k(u^,  u^ )  (surface  constraint). 

u^  =  f(X)  (parametric  separation,  first  variable). 

^2  =  g(^.  a,  b)  (parametric  minima  solution,  second  variable). 

As  =  s(X,  a,  ta)  (parametric  distance  traveled). 

ip  =  Tp  (X ,  a,  b)  (azimuthal  heading,  0-360,  off  a  known  constant/known 
direction) . 

where  X  is  the  parametric  separation  variable,  and  a  and  b  are  constants  of 
integrations  which  (assumably)  uniquely  determine  the  curve. 

Several  surfaces  are  addressed,  but  most  rigorously  treated  is  the  sphere, 
where  the  second  order,  nonlinear  differential  equation  which  governs  the 
extremum  solution  is  solved  for  both  boundary  condition  and  initial  condition 
cases,  as  both  cases  are  of  physical  interest. 


THIS  TECHNICAL  NOTE  IS  THE  WORK  OF  OR.  K.  CARLTON-WIPPERN  ANO  DOES  NOT 
NECESSARILY  REPRESENT  THE  OFFICIAL  VIEW  OF  THE  US  AIR  FORCE. 
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The  nature  o-F  the  oblateness  of  the  earth  is  very  briefly  discussed,  and  a 

first  order,  perturbation  approach  (based  on  e  ,  the  square  of  the  eccentric- 

2 

ity)  set  of  equations  are  developed.  A  comparison  of  the  e  perturbed  solu¬ 
tions  to  those  of  the  unperturbed  spherical  solutions  is  discussed  in  appendix 

!•  The  e  perturbed  solution  set  constitutes  the  mathematical  foundation  for 
Project  Frost  (First  order  formulae  relating  to  oblate  spheriod  trigonometry). 

While  the  purpose  of  this  document  is  to  provide  useable  products,  not  just 
address  theoretical  development,  some  attention  has  been  paid  to  the  existence 
and  uniqueness  of  solutions,  particularly  in  the  spherical  case,  as  there  do 
exist  traps  for  the  unwary. 

The  application  for  the  equations  and  solutions  for  the  problems  addressed  in 
this  document  runs  the  gamut  from  commercial  to  governmental,  particularly 
mapping  agencies,  and  tactical  and  strategic  defense  organizations. 

The  mathematics  developed  herein  are  for  use  only  as  solutions  to  problems  of 
surface  geometry,  curves  of  minimum  distance.  No  attempt  has  been  made  to 
correct  for  physical  phenomena,  such  as  wind,  geopotential,  tides,  currents 
etc.  Corrections  for  heading  and  positioning  are  outside  the  scope  of  this 
study. 

The  format  used  is  necessarily  terse,  and  the  document  assumes  the  reader  to 
be  proficient  in  variational  calculus,  differential  equations  and  linear 
algebra,  as  well  as  spherical  trigonometry  and  in  other  applied 
mat hematics/ physics  curricula. 
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I.  INTRODUCTION 

The  investigation  begins  with  the  general  distance  equation: 
(ds,ds)  =  g^j^dx^dx^ 


or 


s 


gag 


dx 


a 


dt 


dt  / 


dt 


for  a  space  curve  given  by  x^  =  x^(t}. 

The  curves -of  minimum  distance  (or  geodesics]  are  found  by  solving: 


2  X 


ds 


+ 


a  I 
pX  dx  dx 

I  „  ds  ds 

ag 


Q  . 


where 


,  _  JLif  !!2£  ,  ^ 

1 

kind. ' 

For  the  purpose  of  this  investigation,  consider  curvilinear  coordinates  only 

1  2. 

(which  greatly  simplifies  things,  reference  Lass  or  MorsfPand  Feshbach  ) 

2 

also,  the  investigation  primarily  works  with  R  ,  and  assumes  if  u^,  and 
Uj  are  curvilinear  coordinates  there  exists  a  surface  function,  k,  such  that: 


^  ,  called  the  Christoffel  symbols  of  the  second 


Ug  =  k  (u^  ,  ) 


The  distance  formula  becomes. 


^Lass,  Vector  and  Tensor  Analysis,  p.  288. 

^Lass,  op.cit.,  p.  281. 

2 

Morse  and  Feshbach,  Methods  of  Theoretical  Physics,  p.  47. 
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3 

where  h^ ,  h^ ,  and  h^  can  be  called  scaling  factors. 

Seek  parametric  solutions  for  Uj^,  U2  which  will  minimize  s,  since 


U3  =  k(u^ ,  U2  ) 

The  solutions  (for  u^,  U2)  can  be  written  parametrically  of  the  form: 


=  f(X]  , 

□2  =  g  (A  ) 

f  can  be  chosen  arbitrarily  (to  facilitate  integration),  find  g(X}  via  the 

4 

usual  methods  of  variational  calculus. 


“  ■  X  hjl,— 

1=1 


with  u^  =  k(f(X),  u^),  u^  =  f{A),  then  extremum  is  when: 


d 

dX 


86 


3  (U2} 


36  n  •* 

IT -  =  0,  u^ 

3u„  2 


du, 

dX 


letting  U2  =  g. 


L  .  ^  »  L-( ^ 


dX  9  g  \  dX 


dX 


1 

T 


^  ’2 


2  VdX 


h.  •  -Ev 


^  8 

dX 


8g 


dX 


=  0 


(2) 


Rearranging  2a,  taking  a  degree  of  freedom  to  facilitate  integration,  for 
example,  assume  there  exists  a  function,  f  (=  u^),  which  sets 


d 

dX 


n  A  2  dk  8  /dk^\1  1  .2  ^  L.f^\=n 

[6  V^3  *  dX  *  3gVdXyyJ  6  *  ^3  *  dX  8gVdXy 


(2a) 


^Spiegel,  Advanced  Mathematics  for  Engineers  and  Scientists,  p.  288. 
^ Courant-Hilbert,  Methods  of  Mathematical  Physics,  p.  164. 
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then  taking  this  degree  of  freedom  for  f  leaves 


(2b) 


h  g  =  c6 


^4  '1  2 

h^g  =  c 


but 


^2  (  df  Y  f  dk 


(3) 


^  ak  ^  " 

dX  3X  ag  ^ 


4  '2 


2.  2  '2 


0  h„  g  =  c 


2  A  df 


1  \  dX 


+  h; 


ak 

ax 


Y  f  ak  V  '2  „  8k  ak  -■ 
j  *  ( 3?  j  2  *  2  -jj-  .  ^ 


Rearranging  and  collecting  terms: 
-2 


Ag  ^  +  2Bg  +  C  =  Q 

A(X,g)  . 


.  (3a) 


(4) 


„ ,,  ,  2,2  (  ak 

B(X,g)  ^ 


Ui 

9g 


C(X,g)  =  0^ 
Now  calling 


2  /  af 


1  V  8X 


+  h. 


aj< 
3  V  3X 


n(x,g)  =  -  -^<16  ±  p  -  Ac]^/^j> 


And  assuming  n  can  be  uniquely  determined  by  the  boundary,  surface/other  logi¬ 
cal  elimination  conditions,  then 


g  =  n  (X,g)  . 


(5) 


The  geodesic  curves  are  found  by  choosing  f  to  make  equation  3  possible,  then 
solving  the  integral  equation: 


g(X) 


n(X,g)  dX 


(6) 


This  can  be  considered  a  navigation  kernal,  in  that  for  a  surface  defined  by 
k,  the  constraint  condition  for  f,  equation  (2a)  and  the  solution  to  equation 
(6)  gives  us  the  parametric  equations  of  the  path(s)  of  minimum  distance,  sub¬ 
ject  to  the  boundary  conditions  (end  points  of  the  curve  on  the  surface). 
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While  traveling  on  a  curve  of  minimum  distance,  it  is  convenient  to  know  in 
what  direction  to  go  to  remain  on  the  curve.  Using  the  natural  definition  of 

the  tangent  vector  ("=  ^')aids  getting  around  on  the  curve.  If, 


u 

u 

u 


1 

2 

3 


D 


f(X) 

k(X),  and  if  e^^ ,  e^,  are  curvilinear  unit  vectors,  then: 

3  f  9g  3k 

3A  *^2  3A  ^2  ^3  3X  * 


(Compare  this 


result  to  the  contravariant  vector. 


5 


) 


If.  the  surface  on  which  the  navigation  is  taking  place  can  be  written  as 
S(uj,  U2,  U3 )  =0,  and  the  parametric  curve  of  minimum  distance  can  be 

written  also  as  a  surface  ?(u^,  u^,  u^)  =  0  such  that  the  two  surfaces  inter¬ 
sect  in  the  space  curve  of  minimum  distance,  then  the  heading  vector  H  can  be  • 
written  as : 


H  =  Vi;  X  VS  . 

To  demonstrate  this,  let  S,  the  surface  on  which  the  navigation  is  taking 
place,  be  represented  as  locally  (approximately)  planar. 


Similarly,  let  ?,  the  surface  associated  with  the  minimum  distance  constraint, 
be  approximately  planar  (this  is  roughly  equivalent  to  using  affine  approxima¬ 
tions  to  the  actual  surfaces  about  a  small,  differential  area  of  local  inter¬ 
section,  assuming  the  intersection  exists). 


C 

Lass,  op.cit.,  pgs.  270-272. 
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intersection 


The  intersection  (of  the  two  planes)  is  a  line.  Let  N  be  the  surface  normal 
to  S  and  be  the  surface  normal  to  the  direction  of  H  (if  not  i  0)  is 
clearly  x  construction.  In  the  limit  that  the  differential  planar 

surface  approximations  go  to  zero  ^S,  so  that  x 

(gradients  to  be  evaluated  at  a  point  of  interest  on  S). 

So  long  as  s  0,  the  magnitude  of  |H|  is  relatively  unimportant.  It  is  the 
direction  of  which  matters.  Eventually,  |?i|  drops  out  of  the  analysis 
entirely,  it  is  from  the  components  of  ?i  that  heading  angles  are  determined. 
(The  same  argument  is  true  for  6(X).) 


II  EXANPLE^N  CARTESIAN  COORDINATES 


For  this  example  the  surface  is  defined  as  a  plane, 
z  =  a  +  bx  +  cy 

and  search  for  g(X}  such  that  the  integral  is  extrema.  For  example,  let, 
X  =  X 

y  =  g(X} 

z  =  a  +  bX  +  eg  (X )  , 

then 
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2  /  ,  -  \  2 
+  (  b  +  eg  j 


1/2 


dX 


_  ^ 

dX 


3  6 

With  6  =  integrand,  note  that  ^  =  0;  therefore. 


— 0,  and  after  integrating  once: 
dX  \  3  g 


e  = 


/  ^  c  (  b  eg"  ) 


J  +  g"^  +  ( b  +  eg  ) 


(e  is  a  constant  of  integration]. 


Continuing  with: 

[1  +  +  (  b  +  cg')^J  =  [g  +  c  (  b  +  eg  )  . 

This  is  a  quadratic  equation  in  g  ,  of  the  form: 

Ag  +  Bg  +  C  =  0  (A,  B  and  C  are  constants).  . 

Now,  solving  the  equation  for  g  ,  let 

g  =  □,  assuming  D  to  be  uniquely  determinable  from  the  geometry. 


then 


g(A)  =  DX  +  E 
^d  substituting  back 
X  =  X 


y  =  DX  +  E 

z=FX+C  [=a+bX+c(DX+E)] 

3 

This  is  the  parametric  equation  for  a  line  in  R  . 

If  the  line  passes  through  points  ^  ,  y^ ,  ^  and^x^,  y^, 


then, 

y  =  ( vi  ■  Vo )  ^  ^  yo 
^  -  ^0 )  ^  "  ^0 
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is  an  appropriate  expression. 

Letting  a  =  -  x^ ,  b  =  -  y^,  c  = 

then, 

D  =  ai  +  bj  +  ck* 
and  for  distance 

As  (A]  =  (a^  +  AA  . 

The  parametric  equation  can  also  be  written 


or  salving  individually  for  X,  then  setting  components  equal  to  each  other 

X  “  y  -  2  - 

_ 0^  _  ^  ^  _ 0_ 

a  ~  b  ~  c 
The  forms  are  equivalent. 


*□,  a  tangential  vector  (  the  tangential  vector  is  often  referred  to  as 
dr  D  \ 

T  =  -r-,  so  that  T  =  —  ,  this  is  equivalent  to  saying  ^  =  constant  on  the 

iSp 

reference  surface,  where  is  a  heading  angle  measured  on  the  surface, 
referenced  from  some  known  directional  vector. 
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b.  Minimum  distance  solution. 


Consider  the  following  diagram 


To  find  the  minimum  distance  from  p  to  the  curve  g,  let  be  a  plane  (in  par¬ 
ticular,  the  X,  y  plane).  If  f  is  the  solution  to  the  minimum  distance 
problem  on  a  plane,  f  is  a  line  of  the  form  y=ax+borx=ay+b  , 

then  the  minimum  distance  between  p  and  g  lies  along  a  curve  f,  hence  if 
p  =  (  Xq ,  yQ ^  ,  letting  x  be  the  independent  variable,  then, 

(  As  =(x  -  Xg  (  g(x)  -.yg  y 


and  this  is  extremum  when 


ds 

dx 


=  0. 


)‘( 


r(x)  - 


y 

^0  /  dx 


so  that 

(  "  '  =^0  ) 
dx  g(x}  -  y^^  j 


(7) 


(8) 


Solving  for  x  (on  g)  in  equation  (7),  let  the  result  be  x^,  then  find  = 
g(  x^)  jassuming  g(  x^  )  9^  y^j  and,  constructing  a  line  from(  x^,  y^)  to 
(  ^1'  Equation  (8)  indicates  that  the  line  from(^  x^,  y^  )  to 

(  Xi ,  y^  )  intersects  g  at  right  angles.  There  may  be  more  than  one  solution 

to  equation  (7),  in  which  case  the  absolute  (instead  of  relative)  minimum 
would  have  to  be  found,  assuming  one  existed.  Whatever  the  case,  the  solu¬ 
tion,  by  equation  (8),  guarantees  that  the  line  segment  intersects  g  as  a  per¬ 
pendicular. 
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In  that  a  smooth,  well  behavsd  surface  can  be  represented  as  a  limiting  number 
of  planes,  (ie,  a  surface  S  can  be  approximated  by  n,  lim^_^^  differential 

planar  surfaces  or  affine  surface  approximations)  it  seems  plausable  that  if  g 
=  f,  f  a  minima  curve  on  a  surface  S,  then  if  p  is  a  point  on  the  surface  but 
not  on  f,  and  if  f  is  another  minima  curve  which  goes  through  p  and  intersects 
f  at  right  angles,  then  f  is  the  curve  of  minimum  distance  between  p  and  f. 

As  will  be  shown  in  the  spherical  development,  f  need  not  necessarily  be 
unique  for  the  general  case.  Stated  mathematically 

Let  p  e  f 

and  let  f  =  f(A) 

and 

f  =  f  (X)* 

Than  f  minimizes  J~ ds  between  p  and  ^  if  there  exists  X  such  that 
and 

D(^)  •  D(^  )  =  0  ( f ,  f  evaluated  at  X^  j 


Novj  v^/ith  OCf) 


,  db 
+  h  —  e 
2  dX  “2 


hence : 

,2  da  da  ,2  db  db  ,2  dc  dc  _ 

L  dx  *  ^  L  dx  •  dx  ^  L  dx  *  dx  =  ° 


(apriori  uniqueness  not  guaranteed.) 


It  must  be  remembered  that  this  case  is  concerned  with  solutions  to  second 
order  differential  equations  with  boundary  conditions.  Depending  on  the  sur¬ 
face  geometry  and  the  boundary  conditions,  the  minima  solutions  may  not  be 
unique,  as  is  shown  in  the  spherical  section. 

The  Cartesian  example  is  finished  by  working  out  explicity  the  point  to  line 
on  a  plane  example. 

**Compare  this  with  the  definition  of  the  contravariant  vector,  note  also 
ds  =  (□  •  □)^/^dX. 
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y  =  mx  +  b 


(As  =(xo  -  X  +(yo  -  [mx  +  b]  )  ^ 
ds  ^  g  ^  ^  +  ("m](yQ  -  [mx  +  b]  ) 


dx 
hence. 


x„  +  my„  -  mb 

0  0  ^ 
X  =  - - - ,  call  this  paint  x 


1  +  m 


y  =m(x)+b  = 


mXQ  +  m  yo  +  b 


1  +  m 


Slope  of  intersecting  line 


y  -  Yn 


X  -  x^ 


mx^  +  m  y^  +  b 


1  +  m 


0  mxg  -  myg  +  b 


Xg  +  my^  -  mb 


1  +  m 


■m  X  +  my  ■  -  mb 
0  0 


m 


NOTE:  Slope  of  y  =  m,  slope  of  intersecting  line  =  ~  •  OED. 

Ill  EXAMPLES  IN  CYLINDRICAL  COORDINATES 

a.  Application  to  the  cylindrical  surface. 


s  = 


dA  y  VdA  y  VdA 


1/2 


dA 


For  the  first  example  define  the  surface  constraint  to  be  that  of  a  cylinder. 


r  =  a 
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Search  for  9  =  g(A.)  such  that  the  integral  is  minima. 
Letting  z  =  aX,  obtain 


/6dX 


where. 


^ 


2I1/2 


Noting  3^  =  0  , 


then. 


d  /86 


dX  vag 


=  0,  and  integrating  once 


Y,  '2\  -1/2 
g  (  1  +  g  =  C 


'2  2  2  '2 
g  =c+cg  ,g  = 


(1  -  0=^  ) 


2  X  1/2 


g(X)  =  gx  +  0Q  .,0  (X)  =  0^  +  3X 


The  solution  in  parametric  format 


z  =  aX 


9  =  9q  +  gX 


also  z=^(9  -^o)  eliminating  xj . 


1  *  solves  as 


As  =  a  (  1  +  •  AX 


Letting  e^  =  ,6^=1^  ,  e^  =  i^ 
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D  =  ai  +  r3i.Q  =  ai  +  a3ia 

2  y  Z  o 

NOTE:  0  =  Constant,  independent  of  X. 


0  •  i. 


^  is  defined  as  tan  = 


=  g,  ie,  4'  is  a  constant. 


0  •  i 


Similarly,  the  surface  function  is  defined  as  S:  (r  -  a)  = 
distance  function  is  f:  “Z  f  9  “  9g  ^  =0, 


=  0  and  the  minimum 


3 


then. 


H  =  Vf  •  7S 


Vf  =  -1i  + 

z 


Kf)^ 


"0  = 


■^z  3^0 


VS  =  i  hence, 
r 


0  1 


-1  0 


1  ^ 

3  ^z  ^0  • 


This  results  in  tan  ij;  = 


H  •  i. 


H  •  i 


=  3. 
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NOTE:  The  heading  angle  is  easily  obtained  when  either  0  or  H  can  be 
written  in  only  two  coordinate  directions,  that  is  to  say,  given  curvilinear 
u  ,  u  ,  u  ,  if  the  surface  can  be  written  such  that  u,  =  constant,  ip  is  easily 

found  by  dot  product  of  0  or  H  by  the  remaining  unit  vectors  e^,  et  (i,j  t  k). 

^  ^  J 

If  not,  rotate  D  or  H  as  necessary  to  produce  a  meaningful  result  (reference 
IVb  for  additional  data  on  the  heading  angle]. 

The  five  equation  solution  for  the  cylinder  is: 

r(X)  =  a 

z  (A  ]  =  aX 


0(x}  =  9  ±  ex 

0 


As  (X  ]  =  a  (  1  +  0^  )  ^^^X 


(X )  =  tan  ^  0 

b.  Application  to  the  spherical  surface.  Another  example  in  cylindrical 
coordinates  is  as  follows: 


J 


f  6r\^  ^  if 


V  dX  y  V  dX 


+  r 


dX 


1/2 


dX 


The  surface  constraint  is  defined  to  be  that  of  a  sphere. 


2  2  2 
r  +  z  =  a 


,  /  2  2  \  I/2I 

,  [z  =  ±  (  a  -  r  ;  J  . 


9 

Choosing  r  =  - —  ,  which  is  obtained  from  the  following  analysis; 


1  +  A 


From, 


5  = 


dz 

dX 


dr 

dX 


+  r 


2  /  ^ 
dX 


1/2 


dX 


and  since  z  =  f(r)  only 


s  = 


1  + 


dz 

dr 


1/2 


dX 
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2  2  2 
Then  with  r  +  z  =  a 


dz  ^  _  r 
dr  ~  z 


2  2 
a  -  r 


222  d0  ^ 

Letting  r=af  ,r  =af  ,e=g  .dT'S 


,'2  2  '2^  1/2 
T  +  f  g  >  dX  . 


1  -  f 


Define  the  integrand  as  6,  and  by  inspection,  —  =  0,  so,  already  integrating 


86 

once,  — =  c 

9g 


'  r  ^2  '2 

f  g  =  c  - -  +  f  g 

1  -  f 


'2  -4 

■2  2  f  f  „-2  -2 

=  c  - -  +  f  g 

U  -  f  J 


Take  a  degree  of  freedom  by  choosing  f  such  that. 


'2  “4 

f  f 


2  U/2 


(i  -  i^y'^ 


=  dX  ,  let  f  =  sin  u,  df  =  cos  udu 


CSC  udu  =  dX 


-cot  u  =  X 


-  (i  -  f" 


1-f^2l  2  2  1 

=  X^  ^  =  a"  +  1  ,  f^  =  . 

r  1  +  X 
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2  2  2  2  a 

Now,  since  r  =  af,  r  =  a  f  ,  r  =  - then  continuing: 


1  +  X 


2  2 
r 


('*'1  ,  a 

/  —  -  f  ^  ^ 

VdX  /  Vdr  dX 
f 6dX 


(l  +  X^ 


s  =  a 

2 


+  (g" ) 


(  1  +  x^  Cl  +  x^  )  ^  1  +  x^  (l  +  x^  ) 

—  =  0  ,  =  0  (integrating  once). 

9g  9g 


'2  2  1  '2  {  .  , 2  \  -2 
g  =C  +  Cg  5  = 


+  ^  g"  ) 
1  +  X 


1  -  c 


g 


3 


(l-3^X^)^ 

g  =  0Q  ±  cos  ^  3X  . 

Now,  parametrically: 

(assuming  r  >  0) 


r  = 


(i 


9  =  0Q  ±  cos  3X 


z  =  ± 


aX 


(i  ■ 

To  eliminate  X,  start  with 

X  =  3  ^  cos  =  3^  cos  (  9q  “  0 ) 

then  choosing. 


1  +  3“^  cos  (e  -00  ) 


TECHNICAL  NOTE 
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continue  with. 


2  2„-2  2/-  -X  2 

r  +  r  3  cos  ^9  -  0q  j  =  a 

2 

rearranging  and  multiplying  by  5 


Collecting  factors  and  rearranging  to  Cartesian  coordinates,  obtain 

ax  +  by  +  cz  =  0  which  is  the  equation  of  a  plane  confined  to  contain  the 

origin. 

NOTE;  X  =  y  =  z  =  0  always  a  solution. 

c.  Application  to  the  conical  surface.  Another  example  in  cylindrical 
coordinates. 


dX 


1/2 


Define  the  surface  constraint  to-be  that  of  a  right  cone: 


r  =  az  , 
let, 
r  =  aX 


s  =  a 


rr 2  , 2  -2i 

J  t  "  ^  g  J 


1/2 


dX 


as  before 


0  ,  c  , 


9g 


then  rearranging: 

,  4  >2  2  2  1.2  '2 

Xg  =cY  +  cXg 


-/t 


df 


(f"  - 1)*/" 


%  ^  t  -i  sec  ^  f 


(X  =  fc  ,  dX  =  cdf) 
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hence: 
r  =  aX 


6  =  00  +  Y  sec"^  i 


As  =  a 


^  2  2 
^2  .  c  y  _ 


.2  2 
A  -  c 


1/2 


dA 


(A  >  c) . 


IV  EXAMPLES  IN  SPHERICAL  COORDINATES 


a.  Application  to  the  conical  surface. 

(ds  =  (dp)  ^  +  p^  (  d(}))  ^  sin^.j)  ( d9  )  ^ 

The  surface  constraint  is  defined  to  be  that  of  a  right  cone, 
^  =  constant,  which  can  be  defined  as 

(()  =  tan  ^  a 
Letting 
p  =  aX 


and 

.  2  ,  2 
sin  (j)  =  Y 


1  +  a 


continue  with : 

3  =  aj^(l  -+  dA 


17 


1  March  1988 
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2,2/  2,2  2  \  '2 

=  YA(yX  -  c  )g 

^  c  /  2,2  2A  -1/2 

=  ±  Y  X  -  c  ; 


Q  .  1  -1  yX 

=  9„  ±  —  sec  — 

0  Y  c 


(1)  Conic  distance  measurement, 


,^2,  2  2 
d  X  -  c 


dX  =  a 


YXdX 
(y^X^  - 


(2)  Conic  heading  vector  development. 


-  tan  a 


^  sec  Y ( 9 


)  +  p  (choosing  the  positive  direction) 


=  Vf  X  Vs  =  1 


p  sin 


sec  u  tan  u 


(u  =  y(0  -  ej) 


-  sec  u  tan  u  i  +  —  in 

J-  ■  ^  I  p  p  0 

P  sin  <j)  / 
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Now,  looking  on  the  conic  surface: 


cos  '(»  =  1h  I 


tan  = 


ac  sec  u  tan  u 

2  .  X 

P  Sin  <|) 
p  sin  (j) 


(  sin  ij;  = 


ac  sec  u  tan  u 

Now, 

sin  <|)  =  Y  , 

and  from  f: 

ac  sec  Y  (  9  “  9  ) 

- ^ • 

Thus, 

■  tshr  ■  If  (  9  -  *0  )  • 

Now,  to  show  the  equivalence  between  H  and  D,  choose 
p  =  aA. 


(J)  =  tan  ^  a 


e  =  +  i  sec-'  f 


obtain. 


d0 


0  =  ai^  +  p0ix  +  P  sin  (J)  (  -g^  j  ^0  (p  =  aA  ,  sin  (|)  =  y) 


19 


1  March  1908 


TECHNICAL  NOTE 


Now, 


0*1  =  a  =  cos  lO 

P  ' 


Notice  when  the  surface  constraint  can  be  written  =  constant,  the  coor¬ 
dinate  system  can  be  .considered  surface  normal.  In  such  a  system,  when  D  is 
evaluated,  only  two  components  are  present.  These  components  are  along 
orthogonal  topocentric  unit  vectors  on  the  reference  surface.  The  five 
equation  solution  for  the  cone  is: 
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(X ) 


b.  Application  to  the  spherical  surface.  Define  the  surface  to  be  that 
of  a  sphere: 

P  =  a  (constant]  and  search  for  9  =  g(X]  such  that  the  integral  is  a  minima. 


s 


d(t> 
dl 


^  2  .2  /d9^^ 

+  p  sin  ^ 


1/2 


dX 


(1)  Solution  by  parametric  separation  of  variables.  The  differen¬ 
tial  equation  for  spherical  integrating  factor  is  obtained  from, 

^  • 

Let, 


1.  =  f 
9  =  g 

and  obtain: 


Noting, 


.  2 

+  sin 


dX 


ii 

9g 


0  , 


then. 


integrating  once: 

c  =  sin^  fg^ 

V  f  +  sin  fg 


2  2  2  2  ^2  4  ^2 

c^f  +  c  sin  fg  =  sin  fg 


2  f  2 

C  -  +  C 

sin  f 


■  2  „ 
Sin  f 


=  g 


c 


2 


•  2  . 
Sin  f 


+  c 


'2 

_s _ 

•  2  . 
Sin  f 


'2 

g 
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Choose  f  such  that: 

f'  ^  J_ 
sin^  f 

df  ^  ^ 

•  2  ^  ,2 

Sin  f  X 


CSC  f  df  =  ^ 


-cot  ~  ~  Y 


f  =  tan"^  X 


so  that 

(j>  =  tan  ^  X 
Then, 


^  =  0 
dX  ^  ' 


(X  =  tan  (j)) . 


dj) 

dX 


1  +  X^ 


sin 


1  +  X 


0  =  g  (X  ) 

s  =  a  J*  6  dX 
where. 


6  = 
Noting, 

ii 

3g 

d  /  86 

"'Ug' 


1  +  x^  +  X^  (  1  +  X‘ 


.  ,  ,  x2-,l/2 


=  0  , 


=  0 


TECHNICAL  NGTE 
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and  integrating  once  gives 

- —  g  =  c  6  c  =  constant  of  integration 

1  +  X^ 


X^  ^2 


1  +  X' 


g  =  c 


-2 


-2 


1  +  X^  /  (  1  +  X^  ) 


4  '1 


'2  '1.1  f  .  ,2  \ 

=  c  +cX(l+X) 


,4  '1  -2,2  '1  -2.4  '1  '1 

Xg  -cXg  -cXg  =c 


(  1  -  )x^g  ^  -  c  "X^-g  *•  =  c 

-2\,2  -21  -2 
)\  -  C  J  =  C 


[( 


1  -  c 


g  = 


1  c 


[( 


X  1  -  c  X  -  c 


1^  ^1  _  ^'i^ih  ’ 


Setting 


,-2  1  -  C 


-2 


and 


dq  =  ^  dX 


g(q)  = 


-t  r_ 

J  q(  q' 


(  -  1  ) 


1/2 


-1 


g(q)  =  9^  -  q 


g(X}  =  9q  ±  sec  ^  ^  (  with  9q  and  3  constants  of  integration) 


Parametrically,  the  curve  is 
p  =  a" 


<j)  =  tan  ^  X 
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0  =  0g  ±  sec"^  ^  . 

To  gain  further  insight  as  to  the  nature  of  this  curve,  work  backwards  as 
follows : 

0  -  0Q  =  sec"^  ^ 

X  =  g  sec  (0  "  ) 

X  =  - - - IS - — 

-eg  cos  (0  “  8g ) 

X  =  _ ZS _ 

a  cos  0  +  b  sin  0 

where, 

a  =  -eg  ^  cos  0Q 

b  =  eg  ^  sin  0Q 

and  substituting  tan  (j)  for  X 

1 3  n  ^  • 

a  cos  0  +  b  sin  6 

NOTE:  This  equation  relates  <()  and  9  explicity,  ie,  (j)  =  f(9)  or  f  =  9 

Continuing  with : 

a  sin  (f)  cos  9  +  b  sin  (j)  sin  9  =  -c  cos  <|)  • 

Multiplying  this  through  by  the  term  p  and  rearranging  gives: 
a[p  sin  (f)  cos  9]  +  b[p  sin  ^  sin  9]  +  c[p  cos  (j)]  =  0 

Now,  recognizing  the  relationship  between  spherical  and  Cartesian  coordinates : 
X  =  p  sin  <{)  cos  9 
y  =  p  sin  (J)  sin  0 
z  =  p  cos  (() 

This  constraint  equation,  in  Cartesian  coordinates,  becomes: 

ax  +  by  +  cz  =  0  • 
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This  is  a  plane  confined  to  pass  through  the  end  points  (boundary  conditions) 
and  the  sphere  origin.  The  solution  for  set  boundary  conditions  is  unique 
provided  the  end  points  are  noncolinear  with  the  origin. 

NOTE:  X  =  y  =  z  =  0,  always  a  solution  to  the  constraint  plane. 

Since  the  solution  is  confined  to  a  plane  (ie,  linear)  and  the  curve  must  be  a 
GREAT  circle  (ie,  includes  the  system  origin)  the  distance  between  end  points, 
s,  can  be  written: 

s  =  a  A  , 

where  A  is  the  angle  between  the  origin  and  vectors  which  point  to  the  enid 
points.  These  vectors  can  be  called: 

>  r  /N 

p.  =  a  sin  cos  0„i  +  sin  sin  0„j  +  cos  (j)„k 

'^0  L  ^0  0  ^0  0"*  ^0  J 

and, 

Pj^  =  a  sin  <j>j^  cos  0^i  +  sin  (()j^  sin  0^j  +  cos  ())^k 
Note, 

IPo  I  -  |Pi  I  =  a'  , 

and  from  linear  algebra 


therefore: 


cos 

A 

=  sin 

^0 

sin 

1 — 1 

cos 

cos  0 

+ 

1 

sin 

sin  sinij)^  sin  0^ 

+  cos  <() 

0 

cos 

A 

=  cos 

'^0 

cos 

<^1 

+ 

sin 

sin 

j^os 

0Q  cos  0^  +  sin  0Q  sin  ©J 

cos 

A 

=  cos 

^0 

cos 

+ 

sin 

<^0 

sin 

cos 

(%  -  b)  • 

This  is  exactly  what  Naperian  spherical  trigonometry  (law  of  cosines)  tells 
the  reader  to  expect.  The  spherical  law  of  sines  can  be  recovered  via 
geometrical  construction. 
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If  a  boundary  value  problem  is  solved  in  such  a  way  that  a  minima  curve  goes 
through  two  points,  p^  and  p^^,  then  the  minima  solution  is  a  plane  of  the 
form 

N  •  X  =  0 
where, 

N  =  (a  b  c)  , 


N  is  the  planar  normal  vector  (plane  of  the  great  circle)  of  arbitrary 
magnitude.  Then  given  p^ ,  p,  N  can  always  be  constructed  (of  arbitrary 

magnitude)  as  N  =  Pq  x  p^  .  If  N  =  0  when  pj^  =  ±  p^  ,  the  solution  is 
nonunique.  In  this  case,  there  are  infinite  solutions  which  can  be  obtained 
by  finding  N  such  that  either 

-j*.  -V 

p^  •  N  =  0  or,  equivalently 

p^  •  N  =  Q 

(1)  Spherical  detail  on  H,  the  heading  vector.  Define  p^  and  p^ 

as  points  on  a  given  minima  curve.  On  the  sphere,  prerotate  the  coordinate 
system  such  that  Pq  is  in  the  same  plane  as  the  zero  longitude  and  refer  to  p^ 

in  terms  of  ,  A9  )  . 

A  review  of  spherical  coordinates. 
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With  Pq  at  0  =0,  and  using  unit  vectors: 
Pq  =  a  (xi  y j  +  zk  ; 


=  sin  ({iiji  +  cos  (j>Qk 


and  for  arbitrary  p, 
P 


i —  =  sin  (|)^  cos  A9i  +  sin  (|)^  sin  A6j  +  cos  (|)^k 


I 

Now, 

-k. 

P, 


0-  X  -r^ 


iPo  I  iPi 


also. 


1  (-  sin  (j)  sin  A9  cos  <j)  )  +  j  (sin  ((>  cos  A9  cos  (j) 
10  10 


sin  (|)q  cos  <|)j^)  +  k  (sin  (j)^  sin  (j)^  sin 


A9) 


>4. 

Vf 

_  Po 

IVfi 

Po  IPo  i 

since 

Vf  = 

^  (xi  + 

a 

and 

Vf 


|Vf| 


Po 


=  sin  <j)j^i  +  cos  (|)j^k 


Therefore : 

H  =  N  X  Vf* 

N  =  Pq  X  p^ 


*By  definition,  also  refer  to  geometrical  interpretation,  since  N  =  V?,  the 
minima  surface,  (ie,  plane) . 
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sin  (|)j^  cos  A9  cos  (j)^  -  cos  (|)j^  sin  (jig  cos  (t)^ 


H  =  sin  <j),  sin  A0 
r  1 


cos  (|)j  sin  (|)q  -  sin  ij)^  cos  A9  sin  Ijq  cos 


This  vector  now  needs  to  be  rotated  into  topocentric  (SEZ)**  coordinates. 
This  is  done  by: 


cos  <t>Q  0  -sin  cIJq 
0  1  0 


sin  (])„  0  cos  <pf 


calling, 
c  =  cos  (j), 


s  =  sin  <|), 


X  =  sin  <t)j^  cos  A0 


y  =  sin  sin  A0 


z  =  cos  (j). 


3  2  3  2 

xc  -  zsc  -  zs  +  xs  c 


2  2  2  2 
xsc  -  zs  c  +  zcs  -  xc  s 


2  2 
c  (xc  -  zs)  +  s  (xc  -zs. 


**SEZ,  a  local  right-handed  coordinate  system  which  points  south,  east  and  up 
(or  Z). 
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xc  -  zs 

y 

0 


and,  finally,  translating  (using  north  as  a  cardinal  point  and 

substituting] 

/cos  ((ij^  sin  -  sin  (1)^^  cos  (j)^  cos  A9  \ 


NEZ 


sin  <j)j^  sin  A0 


Q 


NOTE:  The  z  component  is  zero,  this  must  always  be  true  in  any  case. 
Now  looking  at  the  .  local  system: 


tan  ^  =  I 

sin  0.,  sin  A0 

tan  ^ 

By  leaving  the  result  as  a  fraction,  ie,  tan  =  •^,  all  quadrants  are  readily 

resolved.  In  an  absolute  sense,  equation  10  could  have  just  as  easily  been 
obtained  via  spherical  trigonometry: 


(9) 

(10] 
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cos  c  =  cos  <j)g  cos  (|)j^  +  sin  ij)g  sin  (t)j^  cos 
cos  (|)j^  =  cos  <i)g  cos  c  +  sin  (|)g  sin  c  cos  ij^ 


sin  _  sin  A9 
^  sin  c 


sin 
continuing: 


sin 
and 

sin  c 


sin  <ti  sin  A9 
:ion  (12) 


c  sin  ip 
from  equation 

cos  Tp  =  — : — fcos  (j>,  -  cos  (j>- 
sin  (})q  L  1  0 

and  substituting  from  equation  (11) 


11 JJ 


(14) 


cos 


sin  c 


1 


sin  c  cos 


sin 

<i>0 

1 

sin 

‘*>0 

1 

sin 

^0 

cos 

'f'l 

cos  ij)^  -  cos  <j)g  ^  cos  (|)p  cos  +  sin  <1)^  sin  ({)^  cos 


A9  ) 


jj:os  (j)^  (  1 


-  COS^  (|)g) 


sin  .|)g  ^  cos  fi)^  sin  <1)^  ) 
sin  <f)g  -  sin  (j)^^  cos  (j)g  cos  A9 


-  cos  <|)g  sin  <{)^  cos  A9  sin  (})g  ^ 

4) ^  cos  <(>g  cos  A6  ) 


sin  (  sin  4> 


(1! 


Dividing  equation  (14)  by  equation  (15)  gives: 


sin  c  sin 
sin  c  cos  t 


=  tan  4'  = 


sin  <j>i  sin  A9 


cos  (|)^  sin  <|)q  -  sin  4)^  cos  4)^  cos  A9 


The  geometrical  interpretation  of  H  =  N  x  Vf  is  determined  from  the  following 
illustrations : 
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In  this  illustration  curve  i  represents  the  intersection  of  plane  and  the 
spherical  surface.  Clearly  H  must  reside  in  therefore  be  normal  to  N 
(ie,  N  •  H  =  0  always). 


From  this  illustration,  2  is  the  planar  surface  normal  to  the  point  p.  H 
must  also  reside  in  this  plane.  Headings  (ie,  north,  south,  east  and  west] 

are  locally  defined  in  this  plane.  H  must  also  be  normal  to  M(M  •  H  =  0 
always) . 

-W  -V 

NOTE:  N  is  the  surface  normal  for  plane  M  is  the  surface  normal  for  plane 

a' ) . 

H  could  be  defined  as  the  intersection  of  and  or: 


H  =  n 
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Since  the  magnitude  of  H  is  not  a  factor,  just  the  direction,  then  the  fact 
that  •  H  =  0  and  =  0  can  be  used.  From  vector  analysis  a  convenient 

way  of  finding  a  vector  %  which  is  normal  to  both  t  and  5  is; 

A  =  C  X  □ 
then  define  fi, 

X  ^  . 

Now  ^  is  obtained  from  and 

^  =  Pq  X  Pi  , 

and  the  surface  normal  vector  is 

=  ^f  , 
where, 

f  =  f(x,  y,  z]  which  defines  the  surface. 

Hence, 

H  =^Pg  X  Pi^  X  ^f  (for  a  sphere). 

^f  to  be  evaluated  at  a  point  in  question. 
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The  space  curve  directional  vector  3  is  found  as  follows: 


For  first  quadrant  direct : 
p  =  a 


=  tan  ^  X 


0  =  9q  -  SBc'^  i 


3(X}  =  Oi  + 


(A  decreasing,  hence  -1  •  dA ) . 
A 


(-l)i,  *  (-1) 

9 


( 1  + 


X 


-  1 


3^ 


Measure  from  true  topographic  north  (ie,  -i,  direction] 


3  •  =  oos  |3  I 

3  •  1q  =  sin  |3  1 


tan  <1)  = 


(i 


1/2 


-  1 
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NOTE:  From  this  cot  =  cos  (|>  tan  (  Oq  “  ®)  • 

l/2\ 


4;  =  tan  ^  f  B  f  ^ 


-  3^ 


For  H  , 

X 

Now,  from  the  surface  constraint: 
P  -  a  =  0  , 

and  for  the  minimum  constraint: 


3  sec(0g  -  9)  -  tan  (|)  =  0  (since  in  the  S  derivation,  a  quadrant  in  which 
0  <  9^  was  used,  thus  making  the  argument  positive). 


= 


IP 


^  sec^  ^  i.  * 


-1 


<{>  p  sin 


(3  sec  A9  tan  A9)  i^  jevaluated  at  p  (j>  ,  A9  ) 


Then , 


sec  <))  _  g  sec  A9  tan  A9 
p  p  sin  (j> 


0 


0 


f,  n--  3  SBC  A0  tan  A9  «  ^  sec  (j)  ? 

H  =  U1  ^  ; - 1  .  ■  ■  ■■  1  -I-  — — • 

p  p  Sin  (j)  (j)  p  y 
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Now,  in  the  local  system 


IH 


cos  ^  =  -cos  ij; 


sin  K  =  sin 


therefore. 


fH 

H 


-cot  ^ 


3  sec  A9  tan  A6 
P  sin  (j) _ 

sec  (|) 

P 


and, 

,  ,  3  sec  A0  , '  .  „ 

cot  4/  =  — r - 7 -  cos  4)  tan  A9 

^  tan 


but. 

3  sec 

A9  =  tan 

<P 

cot  -Ip 

=  cos  ij) 

tan  A9 

OED. 
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Finally,  the  relation  for  could  have  been  derived  directly  from  spherical 
trigonometry. 


given  b,  a  ,  find 
b,  a  implies  a,  c,  d,  t 

ij;  =  IT  -  , 

therefore, 
tan  =  tan  ij) 
with 

and 

c  +  b  =  — 

7 


sin 


sin  c 
sin  s 


cos  d  =  0  =  cos  a  cos  s  +  sin  a  sin  s  cos  i|; 
-COS  a  cos  s  =  sin  a  sin  s  cos  ^ 


-cos 


cos  a  cos  s 
sin  a  sin  s 
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sin  a  ^  tan  c 
cos  a  ”  cos  a 

since 

cos  s  =  sin  a  cos  c 
Continuing  with 


-tan  ^  = 


sin  c 
sin  s 


cos  a  cos  s 
sin  a  sin  s 


sin  c 
cos  s 


tan  a  = 


sin  c 


sin  a  cos  c 


c  +  b 


tan  c 

tan  ^ 

cot  ^ 
but 
a  =  (J) 


cot  b 

cot  b 
cos  a 

cos  a  tan  b 


b  =  A0 

cot  =  cos  (p  tan  A0 


or 


tan  Tp 


NOTE ; 


d  = 


sec  (j) 
tan  A0 

cos  s  = 


cos  a  cos  d  +  sin  a  sin  d  cos  c 


cos  s  =  sin  a  cos  c  for 


for  distance  measurement; 


As 


1/2 


dX 
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As  =  a  r - • 

J  1  ^ 


1 


"4  - 1  *  1 

3 


-il/2 


dX 


As  =  a 


1  .  X^  /  X^ 


1/2 


XdX 


+  1 


1/2 


As  =  a  (  1  +3' 


Substituting: 


1  +  X^  =  u 


2X  dX  =  du 


XdX  . 


X  =  u  -  1 


As  = 


du 


2  J  u(u  - 

Substituting  back  for  u 

(  1  +  3^  )  2  4.-1 

.3  =  a  L_4J_ 


a  =  (3^  +  1  ) 


-(3^  ^  1) 


2 

e  +  1 


1/2 


“22 
1  +  X  -  3 


2 

3  +  1 


1/2 


As  =  a  tan 


-1 


2  2 
X^  -  3 


1  +  3 


1/2  Xi 


0 


The  five  equation  solution  for  the  sphere: 

P  =  a 


=  tan  ^  X 


9=9^1  sec"^  ^ 


As  (X  }  =  a  tan 


-i/x^  -  3"  I  ' 


1  + 


X, 

I  J 

X, 
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t(X)=  tan 


3 


.  1  V'1 

X^  -  / 


From  the  five  equations  for  the  sphere,  combine  the  equations  for  As  and  9, 
eliminating  the  parametric  variable  X. 


9  =  9^±  sec 


X 

3 


As  =  a  tan 


X^  -  3^ 
1  +  3^ 


1/2 


A  g  av 

From  the  spherical  geodesic  (great  circle  route)  -g—  =  the  central  swept 
out  angle,  evaluate  from  the  top  of  the  curve  down  one  side  to  X  =  X  . 


north  pole 


tan 


As 

a 


fLijLi! 

\  1  +  3^ 


1/2 


hence, 

X  =9  sec  A9 


and 
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Note  that  for  a  given  geodesic  3  is  related  to  the  minimum  colatitude  value 
that  geodesic  curve  will  attain.  If  =  $,  then  3  =  tan  $  is  related 

to  the  inclination  of  the  geodesic  curve  (remember,  the  curve  is  coplanar)  by 


hence, 

3  =  cot  i  , 

and 

tan  -  A^i'  )  =  cos  i  tan  (  ^  "  A0  ) 

cot  AtJ;  =  cos  i  cot  A0 
tan  Aij;  •  cos  i  =  tan  A0 

This  is  completely  consistent  with  Napier’s  findings  and  is  sometimes  referred 
to  as  the  spherical  law  of  tangents. 

Although  most  of  the  attention  has  been  northern  hemisphere  oriented,  the 
solution,  by  symmetry,  is  equally  valid  for  southern  hemisphere  trajectories 
■(-<»<  A  <  0].  Singularities  occur  at  X  =  0  (either  north  or  south  pole)  and 
at  X  ->■  “  at  the  equator.  From  a  computational  point  of  view,  equatorial,  near 
equatorial  and  transequatorial  routes  (X  >  “»)  are  inconvenient  to  calculate. 

-1 

Defining  the  variable  x  =  ^  instead  of  X  ->■  «>,  let  x  ^  0.  For  this  case, 
the  formulas  become: 

p  =  a 

Ip  (x)  =  cot  ^  X  =  ^  -  tan  ^  x 


0 (x)  =  ±  COS  ^  (3x) 


As  (x) 


tan 


X 


1  3^  \ 

1  -  3^^  j 


(x) 


1  -  3^x^  / 
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For  geodesics  which  cross  the  equator  note  that  tan  ij;  =  $  at  the  equator.  The 
geodesics  are  assumed  to  be  smooth,  continuous  curves  so  that 'for  a  given  3  in 
the  northern  hemisphere,  the  southern  hemisphere  continuation  of  the  same 
geodesic  will  have  the  save  value  of  S  in  the  absolute  sense,  thus  both  a 
given  geodesic  curve  and  its  slope  are  continuous  across  the  equatorial 
boundary. 

Other  singularities  of  interest  are; 
lim  3^0 

This  becomes  a  transpolar  trajectory; 

7T 

9  =  9q  “  —  on  one  side 


Q  =:  9^  +  on  the  other,  (ie,  0(X]  becomes  a  discontinuous  step  function, 
becoming  undefined  at  X  =  D,  at  either  pole.  Remember  X  >  S  always). 

TT 

lim  3  .  This  becomes  an  equatorial  trajectory;  (j)  e  — .  In  that  A  >  3 

this  becomes  computationally  difficult  to  handle.  The  matrix  solution  removes 
this  singularity. 

lim  A  ->■  3 

(2 

- — 

A^  -  3 

This  can  be  rewritten  as 


i()  (A  }  "  ^  ~ 


*  1  ! 


Thus,  as  A  >-3, 


1);  (A  =  3  )  =  I  . 


For  a  typical  boundary  value  problem,  given  p^,  p^,  find  the  geodesic  curve 
which  contains  Pq  ,  P^.  The  expression  N  =  p^  ^  convenient  to  use.  In 

the  case  of  physical  interest,  where  iti  s  Q,  p^  =  -p^  there  are  an  infinite 

number  of  geodesics  which  will  satisfy  the  condition.  (This  is  not  generally 
true  in  the  slightly  oblate  case,  as  will  be  seen  in  IVb(2). 
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(2)  Existence  and  uniqueness  of  solutions.  One  can  use  the  idea 
of  breaking  up  the  surface  into  small,  differential  grids,  and  numerically 
find  a  curve  between  end  points  through  the  patchwork  of  grids  which  will 
minimize  the  overall  distance;  but,  in  general,  there  is  no  guarantee  that 

A,  ^ 

this  curve  is  unique.  Neither  for  the  curve  f,  nor  f  (as  defined  in  the 
closest  distance  planar  section).  This  is  most  dramatically  seen  on  the 
sphere.  Although  there  are  an  infinite  number  of  cases,  only  two  are  cited. 

(a)  Start  point  north  pole,  end  point  south  pole.  The  minima 

■A. 

curve  exists,  but  is  not  unique.  Take  any  meridian  to  the  south  pole  (f 
exists,  not  unique). 

(b)  Point  p  at  the  north  pole,  the  object  is  to  get  to  the 

equator.  The  intercept  curve  is  any  meridian  to  the  equator  (f  exists  and 
satisfies  all  conditions,  but  is  not  unique). 

Uniqueness  is  highly  surface  and  end  point  dependent.  Solutions  on  the  sphere 
show  a  parallel  of  sorts  with  the  general  solutions  of  a  second  order  dif¬ 
ferential  equation  as  they  should. 

Given  a  point  and  a  slope,  the  minima  curve  is  uniquely  determined.  Given  end 
points,  the  minima  curve  may  not  be  unique  and  there  may  exist  infinite 
solutions. 

Looking  at  the  spherical  case  in  more  detail,  two  cases  (or  subcases  of  par¬ 
ticular  significant  physical  interest)  are  identified. 

Case  1,  given  and  determine  the  minima  curve  parameters. 

Case  2,  given  p^^  and  p ,  determine  the  minima  curve  parameters. 

In  each  case,  assume  p^  is  at  0  =  0.  (If  not,  use  relative  9  coordinates. 
Remember,  the  surface  is  one  of  revolution  and  symmetry  can  be  used.) 

And  in  each  case  it  is  most  convenient  to  work  with  the  solution  planar  normal 
which  in  case  2  is  easily  determined  via  =  p^  x  p^ 

Now,  consider  the  relation  of  with  t  : 


*<1^,  a  topocentrically  defined  heading  angle  (0  <  <  360). 
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NOTE:  sin  3  =  N  •  i,, 

l\ 

Rotate  the  normal  up  to  the  point  (at  this  is  a  y  axis  rotation  through 

7T 

an  angle  that 


which  must  then  be  rotated  (-6)  about  the  z  axis  for  proper  inertial  (I,.],K) 
orientation  (0  is  defined  here  as  the  longitude^of  p  ) . 

-y 

This  equation  relates  the  two  cases,  as  given  N,  which  is  equal  to  Pq  x  P, 
and  forcing^  ]N  I  =  1  ,  then  ip  at  1)^  is  determined  via  the  above  relation, 

and  given  p^ ,  t  then  N  is  obtained  (remember  to  keep  track  of  9,  which  was 
assumed  to  be  zero) . 
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Now  the  angle  3. is  easily  determined  via  sin  3  =  *  i|^  Tor  either  case 

(reference  illustration)*. 


NOTE:  9q  the  spherical/sphere  case  parameter  is  obtained  from 


TT  *"1  9  •  i 

from  which  (depending  on  quadrant)  9  =  —  +  tan  - - 

0  2  ^  ,  I 


node 
By  geometry 


it  X  r5. 


^node  along  the  equator  to  where  the  plane  defined  by  N  intersects  the 

sphere  equatorial  plane. 


The  angle  3  is  related  to  the  spherical/sphere  case  parameter  3  via  3  = 
tan  3. 

Given  a  curve  (9„ ,  3)  and  some  arbitrary  point  p  ,  determine  if  the  closest 

0  5 

point  on  the  curve  to  point  p  can  be  found. 

a 


From, 

tan  ())  =  3  sec  (9  “  9^ ) 


obtain  the  curve  in  the  cartesian  form 
ax  +  bv  +  cz  =  0 


('^  •  X  =  0) 


and  let, 

-V  ^  ^ 

P  =  xi  +  yj+zk 
a  a  -^a  a 


Now,  note  if 


ax  +  by  +  cz  =0 
a  a  a 


then  is  on  the  curve, 
a 


If  this  is  so,  stop  the  process. 


If  not,  continue. 


*Require  |^  |  =  1 . 
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Let  l^l  =  1  ,  =  1* 

construct  %,  S  such  that 

A  •  ■§  =  0 
A  +  ^  ^ 

a 


So: 


A  =  •  ^1^ 

S  =  A  -  A  . 

a 

By  geometry  §  points  to  the  solution  point  with 


N 


^  =  P  -  1^  ^  tl\ti  , 

a  a 
then  let 

|B|  =  1 

and  from  this  S  the  solution  is  extracted. 
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NOTE:  If  P  is  in  the  direction  of  5  =  0,  there  is  no  unique  solution. 

Also  if^  *^  =  0,  ^  =  0,  §  =  the  point  is  on  the  curve  (this  is  an 
a  9 

application  of  Gram-Schmidt  orthogonalization) . 

(3)  Solution  via  matrix  methods. 

(a)  Development.  In  that,  if  ?(t)  is  the  position  vector  of 

the  solution  of  the  spherical  case  (say  the  general  case,  where  P($)  is  a 
vector  function  of  time  measured  from  the  origin)  it  is  easily  shown  that 
for  a  geodesic: 

_1_.  ?(t)  •  ?(t)  =  a^,  a  constant  for  any  time. 

2.  There  exists  a  vector  N  such  that  N  •  P(t)  =  0,  for 

all  t. 

The  two  results  indicate  a  solution  obtainable  in  terms  of  vector  matricies. 

Consider  the  following  diagram. 

z(z') 


Lets  say  ^  is  known  in  the  x,  y,  z  systems  and  |fti|  =  1 
N  •  k  =  cos  g 

-N  X  k  =  9  (which  points  in  the  direction  of  x  ] . 
Now 
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y 


Let  5  have  components  ai  +  bj 
hence 


a 


tan 


(all  quadrants  resolved). 


Ihe  known  vector  contains  the  two  necessary  components  needed  to  exactly 
specify  the  curve  (of  minimum  distance) . 

From  ^(t)  •  ^(t)  =  constant  assume  a  solution  in  the  •  x  plane  of  the  form; 


where,  if  the  minimum  distance  trajectory  is  one  of  constant  angular  velocity, 

j)  =  ujt  +  <5 ,  <5  an  offset  angle,  ij)  the  internal  swept  out  angle.  This  vector 
referenced  to  the  x' ,  y  ,  z  system  (a  rotation  of  6  about  the  x  axis)  is: 

/I  0 

^  (t)  =1  0  cos  0 

\0  sin  0 

and  finally  referenced  to  the  x,  y,  z  system. 
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Let  a,  3  be  the  rotation  angles  which  take  ?  in  the  3 
actual  position  vector. 


a  0  plane,  to  ?  the 


/cos  a  -sin  a  0  \  /  1  0  0  \ 

^  =1  sin  ct  cos  a  0  1  (  0  cos  3  -sin  3  |  f  p 


1  /  \  0  sin  3  cos 


cos  a  -sin  a  cos 


sin  a  sin  3  \  /  cos  ip 


sin  a  cos  ot  cos  3  -sin  3  cos  a  sin 


sin  3 


cos  3 


cos  a  -sin  a  cos 


sin  a  sin 


^  =  sin  a  cos  a  cos  3  -sin  3  cos  a  0 


sin  3 


p  =  [cos  ot  cos  -  sin  a  cos  3  sin  efiji  +  (sin  a  cos  iJj  +  cos  a  cos 
+  (sin  3  sin  t)(< 


3  sin  ii;)j 


=  (sin  a  sin  3)i  -  (sin  3  cos  a}j  +  (cos  3)k  . 

If  ^  is  given  (from  boundary  or  initial  conditions)  a,  3  are  easily  deter¬ 
mined. 

The  objective  now  is  to  show: 


=  X  . 


NOTE:  IP  =  1 


,  |l^  X  i^l  =  1,  since  11^1=1,  |l^|=1,f^*N  =  0. 


In  order  to  keep 


sin  ije  •>  cos 
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and 

cos  >  -  sin  ijJ 

for  the  case  of  taking  unit  derivatives. 
Then, 


dt 


-cos  a  sin  "  sin  a  cos  8  cos  ij; 
cos  a  cos  3  cos  ijj  -  sin  a  sin  ijl 
sin  B  cos  ^ 

for  lj  X 


\ 

/ 


k 

N„ 


\ 

/ 


i(N  P  - 
1  2 


P  N 
1  2 


+  j(PN  -NP)+k(NP  -  PN) 
^  0  2  0  2  0  1  0  1 


With : 

Ng  =  sin  a  sin  8 
Nj^  =  -cos  a  sin  8 
=  cos  8 

Pg  =  cos  a  cos  ^  -  sin  a  cos  8  sin  ij5 
Pj^  =  sin  a  cos  +  cos  a  cos  8  sin  ijl 
P^  =  sin  8  sin  tjJ 
Evaluating  each  component 
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X  p)  •  i  =  -cos  a  sin  3  (sin  3  sin  ip)  -  cos  3  (sin  a  cos  ip 
+  cos  a  cos  3  sin  ip) 

2  -  -  2  .  - 

=  -sin  3  cos  a  sin  <p  -  cos  3  sin  a  cos  ip  -  cos  a  cos  3  sin  ip 

=  -cos  3  sin  a  cos  ip  -  cos  a  sin  ip 

(i\i  X  p)  •  3  =  cos  3  (cos  a  cos  ip  -  sin  a  cos  3  sin  ip)  -  sin  a  sin  3  sin  ip 

-  .  2  -  .  2  .  - 

=  cos  3  cos  a  cos  ip  -  sin  a  cos  g  sin  <p  -  sin  a  sin  3  sin  ip 

=  cos  3  cos  a  cos  ip  -  sin  a  sin  ip 

(i^  X  p)  •  k  =  sin  a  sin  3  (sin  a  cos  ip  +  cos  a  cos  3  sin  ^) 

+  cos  a  sin  3  (cos  a  cos  ip  -  sin  a  cos  3  sin  ip) 

=  sin  a  sin  3  cos  ip  +  sin  a  cos  a  sin  3  cos.  3  sin  ip 

2  -  .  .  •  “ 

+  cos  a  sin  3  cos  ip  -  cos  a  sin  a  sin  3  cos  3  sin  ip 


=  sin  3  cos  >p 


hence. 


dP 

dt 


N  X 


-> 

P 


vt  X  ^s 


(as  expected). 


Since  N  =  and  the  surface  gradient  for  a  sphere  =  P,  what  has  been 
shown  is  that 


d^ 


H  which  is  what  would  be  intuitively  expected. 


dt 

Continuing  with. 


(ip  =  ii)  (t)  ±  <S) 
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Finally,  the  term  6,  t  =  0  needs  to  be  evaluated. 

Lets  say  the  term  ^(0)  is  known 
then. 


cos  a 
^(0)  =  I  sin  a 
0 


-sin  a  0 
cos  a  0 

0  1 


Q 


0 


0  cos  3  -sin  3 
0  sin  3  cos  3 


from  which  6  is  easily  solved. 


NOTE:  For  the  argument  6  may  be  iji  =  03(t)  +  6  or  =  co(t}  -  6. 
The  ±  determination  can  always  be  made  from  the  fact  that, 

P  •  5 


COS  6  = 


F.  I F 


sin  6  = 


F„  X  5 


ip 


Direct  trajectory  can  be  defined  for  o)  and  retrograde  trajectory  for  -oj  as 
usual* 

Now,  look  for  heading  information  for  ?(t). 

Given  P(t)  and  with  known  factors  {h>,  0}  or  {latitude,  longitudef  data 

dl^(t) 
dt 


can  be  found.  Then  with  the  {latitude,  longitude}  data  rotate  this 


dP 


d? 


vector  ^  into  the  topocentric  SEZ  system  (the  Z  component  of  —  mus 
zero  in  this  system)  and  obtain: 
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S 


Let, 

=  sin  fj)  cos  0i  +  sin  (j)  sin  9j  +  cos  (i)k 
and,  noting  (j)  +  L  =  —  (L,  the  latitude  variable) 
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so  that 


cos  <j)  cos  9  cos  ^  sin  9  -sin  <j) 


=  I  -sin  0 

1 


cos  9 


sin  ip  cos  9  sin  (|)  sin  9  cos  ()> 


sin  (j)  cos  9 
sin  <J)  sin  9 


cos  <p 


0  =  k  . 


This  is  exactly  as  it  should  be  since  =  k  as  seen  in  the  local  topocentric 
system.  Letting  M  be  the  topocentric  rotation  matrix, 

=  M 

1  0 


sin  (j)  cos  9 


Pq  =  sin  (j>  sin  9 


cos  (() 


cos  a  cos  'p  -  sin  a  cos  g  sin  ip 


sin  a  cos  ip  +  cos  a  cos  g  sin  ip 


sin  g  sin  ip 


a,g,-p 


(this  relates  <p,  9  to  a,  g,  ip.  ip  is  the  internal  swept  out  angle), 


For  the  heading  vector  development  at  the  point  P,  let  H  =  M  H.  H  is  then 
the  heading  vector  in  the  topocentric  system. 

•>  d'^ 

NOTE:  H  •  k  ^  Q,  since  at  the  surface  point  ri^  is  tangential  to  the  surface. 

5p,  the  heading  angle  as  measured  in  the  local  topocentric  plane  can  be  obtained 
as  usual. 


■ip  =  tan 


^r  •  J  \ 

- - X-  (note  the  negative  sign.) 

H  •  i  / 
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/V 

(i) 


(b)  Summary  of  matrix  solution  for  the  spherical  case. 

Let  P  be  the  position  vector: 

/  sin  <p  cos  6 


^  = 


sin  d)  sin  9 


cos  (j) 

From  the  surface  constraint 


From  the  minimum  distance  constraint  obtain 
P  •  N  =  0  . 


Let 


cos 

sin  )  be  the  vector  measured  in  the  ^  =  0  plane  [central  angle, 
0 


,  measured  from  the  equatorial  node).  Then  9^,  3  initial/boundary  parameters 
can  be  translated  to  a,  y  angle  constraints  such  that 


0  , 
0 


(cos  cc  -sin  a  cos 

sin  a  cos  a  cos 

0  sin  y 

3  determine  ti  and. 


y 

y 


sin  a  sin 
-cos  a  sin 


cos  Y 


/ 


cos  t 
sin  ijJ 
0 


\ 

/ 
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li  =  sin  a  sin  Yi  -  cos  a  sin  Yj  cos  Yk 


ip  can  be  written: 


ifCt)  =  K  +  —  At  (for  constant  velocity  trajectory). 

U  0 


cos  a  -sin  a  cos  Y  sin  a  sin  Y 

R  =  I  sin  a  cos  a  cos  Y  -oos  a  sin  Y 

0  sin  Y  cos  Y 


(?i,  the  directional  vector) 
dt 


'('q  can  be  determined  as 


cos  4' 


-  .  N  -(i^.  X «)  . 


X  I 


.  if  Pn  =  1- 


NOTE:  All  development  was  based  on  the  solution  to  second  order  differential 
equations. 

The  sphere  was  treated  with  both  condition  types,  ie. 

Case  1 : 

□E:  Sphere 

BC:  P(t  =  0)  =  ,  ^(t  =  t)  =  1^^  (boundary  value  problem). 

Case  2: 

□E:  Sphere 

> 

IC:  P(t  =  0)  =  ^  (t  =  0)  =  tan  4'  (initial  condition  problem). 

In  both  cases,  it  was  convenient  to  find  1^,  the  planar  normal. 

Both  cases  have  physical  interest  since  case  1  is  the  classical  navigation 
problem,  case  2  can  be  thought  as:  given  the  position  and  direction  of  an 
enemy  aircraft,  where  is  it  going,  what  will  it  most  likely  attack,  minima 
curve  assumed. 

In  the  second  case,  an  uncertainty  in  ip  (ie,  surface  directional  heading, 
where  ij;  =  ij'  ±  AE)  will  generate  a  solution  envelope. 

The  solution  envelope  is  easily  generated  if  is  rotated  to  the  north  pole. 
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The  solution  obtained  is  a  surface  envelope  (or  sectoral  slice)  which  can  be 
useful  in  determining  probable  targets. 

(4)  Selected  problems. 

(a)  Intersection  of  minima  distance  curves.  In  the  plane,  the 
following  expressions  can  be  written  for  the  two  minima  trajectories: 

y  =  ax  +  b 

y  =  cx  +  d 

The  objective  is  to  find  the  minima  trajectories  when 

y  =  y'  , 

X  =  x' 


hence, 

ax  +  b  =  cx  +  d 
(a  -  c)x  =  d  -  b 
therefore. 


X  =  X 


y  -  y  -  a 


d  -  b 
a  -  c 

'd  -  b 


a  -  c 


(a  *  c) 


+  b  =  c 


d  -  b 
a  -  c 


y  = 


a  -  c 


(ad  -  ab  +  ab  -  cd)  = 
(cd  -  cb  +  da  -  cd)  = 


+  d 

1 

a  -  c 
1 

a  -  c 


(ad  -  cb) 

(ad  -  cb)  (results  are  consistent). 
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NOTE:  If  a  =  c  there  is:  no  solution  if  b  ^  d;  no  unique  solution  if  b  =  ,d. 

Now,  no  unique  solution  implies  a  =  c  and  b  =  d,  ie,  the  trajectories  are  on 
the  same  curve. 

No  solution  implies  y  and  y  are  nonintersecting  parallel  lines. 

On  the  sphere  there  exists  similar  cases. 

Let 


N  •  X  =  0  (one  curve) . 

^  •  X  =  0  (another  curve). 

The  vector  which  points  to  the  intersection  of  the  curves  is  easily  seen  to 
be: 


NOTE:  The  intersection  vector  I  must  reside  in  the  plane  iti,  hence: 

1*^  =  0  ,  but  1  must  also  reside  in  so 

t  •  =0  ,  so  that  f  is  orthogonal  to  both  therefore 

t  «  X  ,  now  since  |t  I  =  a  (if  t  f  0),  t  can  be  set  equal  to  x  ,  then 

let,  I  =  1 .  t  now  points  to  the  latitude/ longitude  of  the  intersection/ 
intersections  (remember  there  are  two,  t  and  -t). 

■[  =  ^  X  Kl  (since  I  must  be  a  vector  in  both  planes). 

^  is  unique  if  and  only  if 

X  it  0  . 


If  X  Kl  =  0,  then  either  or  N  =  -N^  which  is  similar  to  case  2  in 

the  plane,  no  unique  solution.* 

There  is  no  analog  to  planar  case  1  (no  solution)  on  the  sphere. 

3 

(b)  Trajectory/ intercept.  In  Cartesian  R  let  the  invader 

parameters  be. 


and  the  interceptor  parameters  be 


*ThG  trajectories  are  on  the  same  planar/sphere  intersection  curve. 
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^^(t)  =  +  '\t 

where  t/0  are  assumed  constants. 

Then  at  some  time  (t  =  x)  it  is  required  that 
Pq  (x  )  =  Pi  (x  ]  . 


The  object  is  to  solve  for  the  components  of  ^i,  given 
Let, 


P 


t/. 


xO 


xl 


Vo'  ^0 


^1' 


’  V' 


V 

V 


zO 

zl 


and  it  is  assumed 


known 


. 


Now 

at  t  = 

=  T 

X 

o 

+ 

V 

xO 

=  ><i 

+ 

V  ,T 

xl 

^0  " 

V  .T 
yO 

+ 

V  ,T 

^0  * 

"zo" 

=  ^1 

+ 

"zl^ 

isolating 

the 

"l 

terms : 

"xl 

^0  ■ 

^1 

xO 

^0  ■ 

yi 

+  V  „ 
yO 

"zl 

■f( 

^0  ■ 

^1 

^"zl 

and 

1  0 


(16) 
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hence  we  have  a  quadratic  equation  in  t 

)  *  (^0  ”  *  '^vo^  )  *  ( 


^0  ■  *  “zo 


''zo^  ) 


.(x„  -  xj^  .(y„  -  yi  x  ( )  '  x  vjx^  x  2( 


^0  ■  Xl  ^xo^ 


2(^0  ■  ^1  ) 


(vj  -  )t'  .  |P„  -  Pj'  X  2(f„  -  ).  J„x 

/  2  2  \ 

or,  assuming  I  Vj^  >  )  this  can  be  rewritten  as: 


at  -  2bT  -  c  =  0 


p- 1^0 

which  has  solutions: 

T  =  ^  jb±(b^  +  ac) 

Assume  x  >  0,  so  the  solution  is  physically  real,  therefor^ 

*  1*0 

The  components  for  \/^  are  found  via  substituting  t  into  equation  (16). 
Thus,  the  values  of  both  and  x,  the  time  of  intercept  are  obtained. 
Now  if  (  <  ^  )  '  ^  physically  real  solution  can  only  exist 

-'h)  •  '^0  °  ■ 

For  example,  let 
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then 


T  = 


^  I' 


(since  \P^  I  ^  then  for  a  solution  to  exist 


2  2 

(The  case  of  Vq  >  ,  the  discriminant  must  be  real  for  a  solution  to  exist.) 


(c)  Curve  of  pursuit. 

NOTE:  Although  this  section  works  the  curve  of  pursuit  problem,  the  intercep¬ 
tion  problem  on  a  sphere  is  also  addressed  here.  This  is  because  the  solution 
technique  is  the  same  (or  very  similar)  for  both  problems. 


yj) 


y  *(x,  y) 


Now, 


dx  o 


-  y 

-  X 


ie,  let 

■  ^10  *  'V' 

=<1  ■  >*10  * 

and  let 

V  =  V.  cos  8 
ox  0 

V  =  V.  sin  0 
oy  0 

then  set  a  parametric  relation  for  x(t),  y(t). 
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^(t) 

V 

r 

where: 

2  2  2 
r^  =  (x  -  x)  +  (y  -  y) 

V  =  \\J\,  the  velocity  (assumed  constant)  of  the  interceptor. 

The  idea  is  the  instantaneous  slope  (or  direction)  of  the  interceptor  points 
along  the  line  (curve  of  minimum  distance)  which  includes  both  interceptor  and 
target  instantaneous  positions.  As  with  most  differential  equations  of  this 

type  exact  solutions  depend  on  [x(t),  y(t)l  the  target  instantaneous  position 
and  the  initial  conditions. 

On  the  sphere,  refers  to  the  target  and  to  the  interceptor.  As  of  the 
target  changes,  iterate  on  ?  and  iti  of  the  interceptor.  It  is  convenient 

(for  N  =  constant  target  on  a  minima  course)  to  use  the  matrix  concepts 
developed  earlier.  Treat  this  case  analgous  to  target 

y  =  y„  *  vy 

s .  .  vt  . 


The  instantaneous  N  is  always  P  x  P, 
divided  up. 


This  is  used  as  the  motion  in  At  is 


Say  for  example  the  given  terms  are: 
calculate: 


K’  P'  (PMVi).  N,  =  Pg  X  Pg  ,  th 


en 


P^  =  M 


cos  (  6  +  I—’  At  j  \ 

a  J\ 


sin  (  6  +  At 
a 


/ 


[remember  M  =  M(N) ;  M,  related  to  the  a,  0  matrix. 


reference  IVb (3) (a)2] . 


63 


1  March  1988 


TECHNICAL  NOTE 


COS  ^5  ±  ^  At 


=  M  sin  (6  ±  L*  ]  At 


M  =  M(NJ 


then  recalculate  N, 


Since  N  is  constant,  M  is  constant  and 


p  =  ri 

n 


COS 

^6  +  n 

|v| 

a 

At 

sin  ' 

f  6  -I-  n 

111 

At 

Now,  since  N  changes  -For  each  At  cycle,  M(N)  changes  and  for  a  given  cycle 
(call  it  M3 


^  X 

m  m  m 


M  =  M(^  ) 
m  m 


CCS  [6  ±  '  At 

V  m  a 


p  .  =  M  sin  (6  ±  I— I  At 

m+1  m  ■  V  m  a 


and  continue  the  cycle. 

As  in  the  planar  case,  the  spherical  case  is  easily  generalized  to  Pq  = 
for  ?’  a  general  curve. 
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NOTE:  The  ^  terms  are  ?  (6  ±  At  )  .  The  ±  can  be  resolved  by  requiring 

^  ^  ^  ->•  ± 

P ,•  P,  >  P„  •  P„  or  P  .  •  P  4  >  P  *  P  ,30  the  interceptor  always  points 

110  0  m+1  m+1  mm  j  ^ 

towards  the  target. 

(Assuming  \\l\  >  |V|). 

(d]  The  spherical  interceptor  problem. 

Given  (and  )  and  (and  find  the  intercept  path. 

Let  be  the  point  of  interception. 

For  the  two  vectors  %  and  S  and  a  unitary  matrix  M 

^  ^  •  m5  (since  only  the  relative  orientation  is  changed). 

Now  let 


^  not  proportional  to  x  has  been  ir 


implicity  assumed.  The 


assumption  should  be  apriori  checked 


by 


i^O  X  ?  J 


* 
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Hence 

(cos 

sin 

0 

cos  6 

M  =  i  sin  6 


(6  known). 


Continue  with 

/  cos  (6  +  ojt)  \ 


M"^Pq(t)  = 


=  I  b' 


sin  (6  +  a)T ) 

\  0 

\ 

(fixed) *. 


(-■w 


Then  find  t  such  that; 


V, 


(t)  •  M  ^Pq  (□)  =  cos  T 


and 


V. 


(target  equation) 


•  ri  ^Pn(T)  =  cos  —  T  (interceptor  equation). 
10  9 


N  ^Pq  (0)  is  defined  as 


I  cos  6  ^ 
sin  S 


,  so  the  first  equation  is  always  true. 


\o  / 

The  second  equation  becomes 


"Plj  the  target  vector  P^,  the  interceptor  vector. 
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a  cos  ^6  ^  ^  ^  ^  ^  ^  \  a^  ^  * 

Now  since  (for  physical  solutions)  x  >  0,  the  x  =  f(x)  type  iteration  solution 
can  be  used: 


T 


-1 

COS 


cos 


+  b  sin  (  6 


starting  with  x  =  0  and  continuing  until  some  convergence  criteria  is  met. 
Then  with  x  solved. 


With  the  point  of  interception  found  known  )the  interceptor  can  be 
navigated  to  the  target  on  the  x  plane. 

(5)  Directional  range  and  location  problem. 

Given  ,  s,  ^  find  (eg,  let  be  a  radio  beacon  at  a  known  location. 

Given  a  radio  receiver  at  v^yhich  can  determine  it’s  distance,  s,  from 
and  the  direction  (i|;)  to  ,  where  is  the  receiver  located,  detsnoine  P^.] 


N 


has  coordinates  (j)|j ,  9^,  and  has  coordinates  (f)^,  6^;  <1)^,  9^,  »j^  and  s  are 
known . 
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sin  \p 
sin  (p^ 


sin  A9 


so  that 


Q  a  •  -1  /  .  sin  ij; 

0.  =  9n  +  sin  (  sin  s  „■  -r 

10  V  sin  <p. 


determined) 


NOTE:  Quadrant  relationship  to  A9  positive. 
Also : 

cos  s  =  cos  4>  cos  <t>  sin  <))  sin  <j)  cos  A0 
0  1  0  1 

cos  IJq  =  cos  cos  s  +  sin  sin  s  cos  4> 


sin  (j) 


cos  S  -  cos  (j)Q  cos 

sin  (})n  cos  A0 


,  X  f  sin  s  cos 

cos  =  cos  S  cos  9,  +  (  — : - - - ^ 

^0  sin  4)q  cos  A9 


^  ^  cos  S  -  cos  cos  (|)^  ) 


/  sin  s  cos  9 

cos  9.  “  cos  s  (  — : — T - rs" 

0  I  sin  9n  cos  A9 


cos  (j>,  (  cos  s 


sin  s  cos  9  cos 
sin  9n  cos  AS 


cos  9 


cos  <j)g  sin  9g  cos  A9  -  cos  s  sin  s  cos  t 
1  cos  s  sin  9n  cos  A9  -  cos  sin  s  cos  t 


^  9  j  determined  ) 


cos  9 


sin  29q  cos  A0  -  sin  2s  cos  9 
1  2  cos  s  sin  9n  cos  A9  -  cot  9n  tan  s  cos  9 


NOTE:  On  an  oblate  planet,  if  it  can  be  assumed  the  locator  beacon  travels 
the  minimum  distance  route,  finding  ,  given  9.  s  is  an  iterative 

process.  Start  with  the  sphere  relations  and  iterate  using  oblate  geodesic 
curve  relations.  Solution  is  unique  as  long  as  is  not  a  pole  {north  or 
south ) . 

cos  s  sin  9g  cos  A9  +  cos  9g  sic  s  cos  9 

but  if  9g  =0  then  cos  9^^  =  cos  s,  A9  undeterminable. 


(6)  Given  range,  calculate  location  problem. 
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(a)  On  a  plane. 

Given  |sj,  [s^l  locate  P^. 


Let 


P.  = 
1 


X. 

1 


and  assume  unique,  so  that 


2  2 


(  -  ^o)  "  (Vl  ■  Vo) 

2  2  2 

(  ^2  "  ^o)  ^  (  ^2  *"  ^0  }  ^  ^2  •  equations  in  two  unknowns  (note, 

there  may  be  two  solutions,  one  solution  or  no  solutions). 

Now 


Xn  =  ± 

0  2 


[=2  -  (y2  -  ^o)  I 


211/2 


so  that 

=^2  ■  -  K  -(^2  -  ^0 

and  solve  for  ^then  for  ^ 
(b)  On  a  sphere. 


s  201/2  ,  02  /  .21 

)  J  n^i  "  (^1  “  ^o)  J 


21  1/2 


69 


1  March  1988 


TECHNICAL  NOTE 


Let 


P  •  P  =  known  =  d 
0  1  1 

^  •  1^  .  =  known  =  d 

0  2  2 

1^  I  =  known  =  1 
'  0  ' 

hence,  assuming  |  =  1 

P  ,  ^  given,  solve  for  ? 
12^  0 


(x,  y,  z  easily  solved  for). 

(7)  Triangulation  problem.  Given  and  find  f,  the 

point  of  interception. 
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(a)  On  a  plane. 


Trivial : 


I 


Let 

1 

—  =  tan  . 

a .  1 


b .  =  y .  - 
1 


a .  X. 
1  1 


where 


hence 

y  =  ^  ( V) 


y  =  a  X  +  b,  f  P  I  )  . 
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agx  +  bo  =  a^x  + 
and 


bi  -  bo 

X  =  -  (  provided  a^  *  a,  ^ 

ao  -  01  ^  0  IJ 


y  =  aQX  +  bo  =  a^^x  +  bj^  (either  relation) 


Hence 


I  solved  for: 


(b)  On  a  sphere. 
Also  trivial,  look: 


As  previously  done  [reference  IVb(2)(b)]  generates  and  f  is  gener¬ 

ated  as  : 


1  .  X 


then  setting  |t  |  =  1  obtain  {<[),  0}  directly.  Always  make  sure  the  quadrant 


is  correct  as  x  points  to  ±  t. 


c.  Application  to  the  slightly  oblate  spheroid.  First  order  perturbation 
techniques  are  used  throughout. 

(1)  Surface  modeling.  The  zero  order  approximation  is  that  of  a 

sphere. 

2  2  2  2 
X  +y  +z  =a  orp=a 
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2  2  2 

Working  with  the  cross  section:  x  +  y  =  a  . 

For  the  first  order  planetary  shape  approximation,  the  assumption  is  made  that 
the  surface  normal  points  roughly  in  the  direction  of  the  acceleration  vector. 

(a)  Nonrotating  surface. 

=  V  (  x^  +  y^  -  )  , 

let 


Sp  /  ^ 

M  .  ^  (xi 


yj  )  »  R  a  scaler  constant. 


Also, 


tan  (j) 


dx 

dy 


or,  the  y  acceleration  component  divided  by  the  x  component. 

2  (b)  For  a  rotating  surface.  The  x  component  is  reduced  by  an 

amount  -co  x  due  to  centrifugal  acceleration.  If  the  resulting  surface  is 
still  to  remain  approximately  normal  to  the  acceleration  vector,  then 
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yielding 

/  2^2  2  ,2/2  co^R 

(  1  -  e  ;  X  -  y  =  b  I  e  =  — 

which  is  of  the  form  of  an  oblate  spheroid. 

Now,  for  a  solid  body 

0  <  g  <  gp, 

since  g  =  g(r]  (assume  symmetric,  isotropic  body). 

g(r)  =  p r  P  =  average  density  (assume  symmetric,  isotropic  body]. 

As  a  first  guess  (and  trying  to  strike  a  compromise  between  a  solid  body  and  a 

/I 

thin  surface)  let  g(r)  =  gp^*  the  arithmetic  average  of  g(r)  on  r,  from 
0  <  r  <  R,  R  the  approximate  radius  of  the  spheroid  volume. 

Then  substituting  back: 

2  30)^ 

<  e  >  a  - . 

2irGp 

For  the  earth 

2  2 
w  2  5.3175*10  (radians/second) 

G  =  6.673*10"^^  N  M^/kg^ 

p  =  5.522  gcm"^ 
yielding 

<  e^  >  2  0.00689  . 

The  plumb  line  is  deflected  as  follows: 

*The  same  value  is  obtained  via  <  gp  > 
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and 

sin  £  ^  sin  (j) 

cos  (j) 


hence 


^  sin  2<p  — 

2  Sr 


0)  R  sin  2(t> 


2g, 


R 


2  7 

Por  an  altsrnative  derivation  of  <  e  >  ,  dimensional  analysis  can  be  used. 
From  the  Buckingham  pi  theorum  choose  the  following  pi: 

77^  =  e,  eccentricity  (a  pure  number). 

7T^  =0),  angular/rotation  velocity. 

TT^  =  p ,  density. 

TT^  =  G,  gravitational  constant  (although  a  constant,  VA/atch!!). 


Start  as  follows; 


Fowles,  Analytical  Mechanics,  P.  130. 

7 

Hughes,  Schaum’s  Outline,  Fluid  Mechanics,  Chapter  5. 
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Assume  w  in  the  numerator  (this  seems  feasable  as  tjf,  et)  then  going  back  to 
the  pi  theorum: 


NOTE:  This  also  assumes  n  >  0. 

1 

If  n  H  is  chosen,  then  the  value  for 

Wstical  "(It)"’' = 

''actual  *  '  ^a  "  ' 

Thus  giving  a  percent  differential  for 

_  3 

(2)  Oblate  spheroid.  According  to  19/5  CRC  for  the  earth, 
a  =  6378.388  km 
b  =  6356.912  km 
yielding 

=  0.00672  . 

NOTE:  These  sections  are  not  meant  to  serve  as  a  rigorous  derivation  for  e, 
the  planetary  eccentricity  value,  but  more  as  a  heuristic  demonstration  on  how 
to  arrive  at  an  approximate  value. 

Consider  the  surface  model  to  be  a  slightly  oblate  spheroid  with  eccentricity 

e  <  1,  e  <<  1.  In  the  analysis  consider  only  terms  up  to  and  including  e  . 
All  other  terms  shall  be  discarded  as  being  vanishingly  small.  Continue  with: 

2  2 

—  +  —  =  1  (letting  r  =  p  sin  ((>  and  z  =  P  cos  <(>)• 

2  2 

a  b 

.  2  2 
Sin  (p  cos  <p  ^  J _ 

2  ,2  2 

a  bp 

^Consolidated  Rubber  Company  (CRC  Handbook  of  Mathematical  Tables,  1975). 


k  is 


>  =  0.083G1  )  . 


)  . 


k,  100 


^t  ■  ^ 


=  1.22%!! 
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Z  Z2/  Z  .2^  2  2^\ 

P  =  a  b  b  sin  ■  <}>  +  a  cos  <j>  ) 

2  ^2  /  b^  .  2  ^  2  ^  \  -1 

P  “  b  I  —  sin  <(»  +  cos  <|)  I 


1  •  2 

1  -1  / 

Sin  (j) 

\  ' 

,  2 

1  1 

.  2  ' 

1 

\-l/2 

'  \ 

1  -  e 

sin  <|) 

) 

n  4-  .  b  2 

let  1 - j  =  e 


Now,  assume  e  <<  1. 

/  0^  2 
p  =  b  M  +  —  sin  (j) 


and 


d(j) 


order  of  e 


set 


dp 


~  0  ^si 


dX 

Continue  'with 


since  ^  ^  ^ 

dX  d(|)  dX 


S  = 


b  I  1  +  -I—  sin^ 


d(j) 

dX 


+  sin 


dA 


1/2 


dA 


As  in  the  case  of  the  sphere  let 


(j)  =  tan  ^  A 


d(|) 

dA 


1  +  A^ 


.  2 

Sin 


so 


1  +  A^ 


s  -  b  1 


2,  2 
e  A 


1/2 


2  (  1  +  A^)  /  \  (l  -f  A^  )  ^  1  A^ 


dA  ( c: 

\  ^ 


Since 


"  dA ; 
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ii=  0  ■ 

9g  U  , 


let 
36 

so  that: 


=  c 


1  + 


2.2 

e  X 


2(1+  X^)_ 


1  +  X 


X2 


=  c 


1»X' 


letting 


c  =  c  1  + 


2,2  \ -1 

e  X  \ 


(1  +  A^) 


p'2  2  /  e^X^ 

c  s  c  1  - 


(1  .x^); 


t  'I-  A 


(  1  +  X^'  )  ^  1  +  x' 


2  N  -2 


,  ‘+  '2  '2  '2.1  (  .  .2\ 

X  g  =C  +cX(l+Xj 

,4-2.2  2  2  X^  2.2  r  ^  ,  2  ] 

X  g  =  c  -  c  e  - -  +  cX  \1+X  / 

1  +  X^ 


2^-2  2  2,4  -2 

g  -  c  e  X  g 


t1  +  c^e^)x‘'g'^  =  c^f'1  -  V  +  X^)g^^ 


1  +  X‘ 


2,  2 
s  X 


X^jx^CI  +  c^e^)  -  c^[1  +  x^j]g"^  =  ^  1  -  — 

/  2  2 

,202  ,222  2  2,  21  -2  2/,  eX 

X  X  +Xce  -  C  -  cXjg  =c|1 - ; 


\  1  +  X' 

/  2  2 

,202,,  22  2,  21  -2  2/,  eX 

X  A(1+ce  -  c)-c  g  =c1 - - 

\  1  +  X^ 
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'2  I 
g  =  1 


2  2 
e  A 


P  c 


\  1  +  X 

(  ^  n2  2  2  \ 

V1"C  +ec  ) 


_  ni/2 


g  =(±)  1  - 


2  2 
e  X 


2(1  .  A^) 


(again  expanding  on  e  ) 


NOTE:  The  second  (e  )  term  acts  like  a  perturbing  term  which 
from  the  great  circle  route. 


Now,  assuming  the  positive  solution: 

■I 

g  = 


e 

1/2  ■  2 


A  _  1 

.2 


TTl 


g(X  )  -  ssc  ^  J  ^  qCX) 


(1  . 

[q(A)  a  perturbing  term). 


q(A) 


■/,7 


AdA 


+  A 


2)  /  ,2 


1/2  • 


-  1 


_et 


y  =  1  +  A' 


du  =  2AdA 


A^  =  y  -  1 

4=1^  (.-1) 

6  3 


q  (A ) 


f- 

2  J  1 


du 


y  (y  -  a) 


1/2 


q(A)  =1 


q(A)  = 


2  -  tan-^^^ 


(a) 


1/2 


(  a  =  1  +  3^  ) 
l/2~ 


(1  +  32)1/2 


tan 


-1  /  X2  -  3^  \  1/2 

1  +  32 


TECHNICAL  NOTE 


always  deti^acts 
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g(X)  =  0Q  ±  |sec'^(|y 


tan 


-1 


NOTE:  g  = 


X  _  -  1 
\ 


/  2  ,2 


So  for  the  first  three  equations: 


P  =  b  1  + 


2  ,2 

e  X 


2  1  .  X^ 


X^  -  3^ 


2  .3^)1/2  \  .  3 


(X  >  b). 


,-1 


<j)  =  tan  1  X 

-1  X  8 

0  ■  2  (TTTT^ 

For  the  sphere  the  nodal  crossings: 


A  A  -I-  1  ■  1  X  B 

0  =  6^  ±  ^  sec  ^  - 


tan 


-I  /  X^  - 

M 


(X  >  B). 


9  ri ■  X  1  =  0n  ±  T  • 
dim  X  ->■  “>)  0  2 


For  the  slightly  perturbed  spheroid 
®  dim  X  >  "  ^0  ■  2  (  ^  ~  2 


(  1  H-  b") 


r 


The  new  curve  is  slightly  crunched  up,  presumably  to  avoid  the  equatorial 
budge. 

The  distance  equation  becomes: 

2 


S  H  b  M1  ^1- 


1 


1  +  x^ 


(l  ^  X^)^ 


xHx^Cx^  -  B^) 


1  -  e' 


2  X' 


1  +  X^ 


l/2 


dX 


s  =  b 


s  =  b 


1 


1  +  X 
1 


/  2,2 

-  1  -H  ^ _ ^ _ 

2  2  .  .  ,2 


d. 


1  +  X 


^)b^ 


1  +  X*-  /  L  X^  -  B^ 


1  -  e 


2  X 


1  +  X' 


1/2 


riX 


1  + 


2  X  2 

e  X 


1  +  X^  (x^  -  B^  )  \  "  1  +  x‘ 


ITz  „2  ^2  -2,2  2,2q21i/2 

X  -B  +B  +  BX  -  eXBI 


dX 
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Combining  terms  yields: 


-  3^)^/n 


(  1  +  3^  )  1 


(l  +3^) 


{X  >  3) 


2\  2  N  1/2 


b  =  b(l+(l-e^)3^) 


(  1  +3^) 


2  ^  1/2 


=  b  1  -  e 


2  3 


2  \  1/2 


1  +  3 


SO  that  (up  to  e  terms) 


As  =  b  1  “ 


^  1+3^ 


.2  ^  1/2 


1  +  —  - 

2 


4(1+3^)/  \  1  +  3 


(a^  -  3^^^^ 


(1^3^) 


This  can  be  reduced  to: 


As  =  b  <  1  + 


■i/x^  -  3" 


4(1+3^)/  \  1  +  3^  /  4(1+3 


,2  U/2 


1  +  A' 


Now  turn  to  the  task  of  finding  3(A)  in  order  to  obtain  instantaneous  heading 
information. 


With , 

p  (A  )  =  b  /  1  +  I - j 

\  ^  1  +  A^ 


(j)  (A  )  =  tan  A 


0  (A )  =01^  +  sec  ^  J  -  ^ 


-llx^  -  3M^/" 


(1  +  3^  ) \  1  +  3 
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dd)  „  ■  a.  d9 

Hr  \  31 


Eventually  the  ratio  of  two  components  must  be  evaluated.  The  vector  l5(X)  is 

2 

divided  through  by  p  and  expanded  (up  to  e  only] : 


2 

Since  terms  up  to  e  must  be  retained,  note  that  the  i^  term  does  not  vanish. 
A  look  at  the  geometry  shows  immediately  what  must  be  done  next. 


NOTE;  i^  does  not  point  along  the  surface,  but  into  the  surface. 
Next,  a  vector  (X)  is  needed,  this  vector  is  normally  measured 


face  (ie,  S(X)  •  s  Q).  What  is  needed  is  a  rotation  about  the  Iq  axis, 
the  angle  is  determined  by  setting  5  •  ip  =  0,  so; 
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First  set: 


TECHNICAL  NOTE 


0  =  -sin  5 


1  +  X' 


+  cos  5 


e^X 


(l  +  X^  ) 


tan  5  = 


e^X 


1  +  X' 


- 1  2 

(0  <  X  <  »  ,  ?  <  tan  (e  )  ,5  bounded  and  small]. 


Now,  since  5  is  very  small,  let  cos  C  ~  1,  sin  K  ~  tan  5 


e^X 


1  +  X 


5'(x]  = 


e^X 


e^X 


1  +  X^  1  +  X^  \(l  +  X^  ) 

U-i-^Vu 


*  (i  '(7^77^7^ 


2  1  +  X 


2  re 


5"(x] 


1 


2  > 


1  ^X-  ^  (1  .  -  3^) 

4 


_  I  _  e^_xi_V- 

2  \  1/2  /  ,  2  „2  \  1/2  2  ,  ,2r0 


3 


1  +  X 


(eliminating  e  terms). 

It  would  appear  from  this  analysis  that  justification  exists  for  keeping  the 
old  ^(X)  and  just  dropping  the  heading  term. 

Now  consider  the  heading  angle  As  usual,  the  work  is  performed'  in  quadrant 

I. 


tan  4) 


2  1 .  (TTTTy^ 


84 


pole 
A0  =  ir 


For  equator  to  equator  trajectories  where  A6  <  rr  ^1  -  ^Jthe  spherical 
geodesic  is  the  same  as  the  oblate  geodesic. 
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For  equator  to  equator  trajectories  where  tt  >  A0  >  tt  ^1  -  ^  j  ,  the  spherical 

geodesic  remains  along  the  equator,  whereas  the  oblate  geodesic  takes  a  route 
similar  to  that  of  the  dotted  line  (actually  there  are  two  geodesic  routes 
(north  and  south)  by  surface  symmetry  about  the  equator]. 

For  equator  to  equator  trajectories  with  A0  =  it  there  are  an  inifinite  number 
of  spherical  geodesics,  there  are  only  two  oblate  geodesic  routes  (one  over 

each  pole).  In  general,  when  =  -p^  ,  for  the  spherical  case  there  are 

infinite  geodesics,  on  the  oblate  surface  there  may  be  only  two  geodesics 
(over  the  poles). 

Thus,  another  effect  of  oblateness  is  to  remove  the  singularity  in  Lhe 
of  meodesics  which  can  occur  (exception:  for  the  oblate  case  there  is  at 
least  one  pair  of  points  which  produce  infinite  geodesics,  north  pole  to  south 

pole) . 

•  (4)  Comparison  of  FROST  to  integrated  great  circle  routs.  Spherical 
minima  routing  gave  a  relation  between  <(.  and  0,  such  that  the  route  was 
confined  to  a  plane.  If  this  route  is  integrated  on  the  slightly  obiate 
spheroid,  the  following  equations  are  obtained: 


p  (X )  =  b  1  +  — 


1  +  X' 


T  (X )  =  tan  ^  A 


-1  X 

(X)  =  0g  ±  sec  j  . 


Then, 


and, 

*Up  to  e^  only. 
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4'*  =  tan  ^  B 


-  3^/ 


Referring  to  the  geodesy  example  in  Appendix  I  (New  York-Leningrad) ,  the 
integrated  great  circle  parameters  are: 


3  =  0.46113733  . 

X  =  1.130294  . 

a 

X^  =  0.577350  . 

The  appropriate  level  of  significance  is  ~  e^ . 

The  FROST  distance  measurement  gave  AS  =  6735.37  km. 

The  integrated  great  circle  measurement  gives  AS  =  6735.37  km. 

This  result  is  not  surprising  if  the  difference  between  the  two  routes  is 
considered: 


In  both  cases  the  value  of  3  gives  the  minimum  colatitude  the  route  will  take. 
For  FROST  -  24.7379°.  For  integrated  great  circle  -  24.7562°. 

The  difference  between  the  two  comes  out  to  be  approximately  1.1  nautical 
miles. 


*Up  to  e^ 


only. 
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Does  this  result  challenge  the  validity  of  FROST?  No,  rather  it  certifies  it. 
This  particular  example  has  starting  and  ending  points  far  north  of  the 
equatorial  bulge  (the  whole  route  is  far  from  the  equator).  It  is  in  cases 
like  this,  where  the  values  of  X  and  g  are  relatively  small  that  one  would 
expect  differences  between  the  two  routes  to  be  vanishingly  small.  Dn  the 
other  hand,  routes  predominantly  olose  to  the  equator  should  show  greater 
distance  and  routing  deviation  (consider  two  points  on  the  equator  179.9° 
apart  in  longitude) . 


K.  Carlton-Wippern ,  PhD. 
Mathematician 
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APPENDIX  I 


10.0  PROJECT  FROST 

10.1  The  need  for  geodesic  and  surface  navigation  algorithms 

There  exists  a  continuing  military  requirement  to  review  and  upgrade  (as 
needed)  the  algorithms  used  to  perform  geodesic  measurement  and  perform  sur¬ 
face  navigation.  These  algorithms  affect  the  manner  in  which: 

a.  Mapping  agencies  measure  distances  between  points  (along  a  geodesic 
or  curve  of  minimum  distance) . 

b.  Long  range  aircraft  (LRA)  navigate  on  the  reference  surface. 

c.  Long  range  naval  operations  (surface/submarine/air)  perform  surface 
navigation. 

d.  Long  range  cruise  missile  surface  guidance  is  performed  during  iner¬ 
tial  flight  phase. 

e.  Strategic  defense  initiative  surface  interceptor  problems  can  be 
addressed. 


(1)  Airborne  platform  positioning  and  navigation. 

(2)  Reentry  intercept. 

(3)  Interception  of  compressed  trajectory  vehicles  via  kinetic  kill 

vehicles . 

f.  Other  (NASA  type  applications). 

10.2  Current  methodology 

Use  spherical  trigonometry,  with  a  compensated  equal-area  radius, 
a.  Basic  spherical  trigonometry  relations. 
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(1)  Law  of  tangants. 
tan  AiJ;  •  cos  a  =  tan  A0 


(2)  Law  of  cosines.  Law  of  cosines  generally  used  to  measure  out 
distances  (can  also  use  tangents  relation  with  proper  parameters). 

cos  Aij;  =  cos  cos  (|),  +  sin  (j)  sin  cj).  cos  A© 
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a  =  b  I  1  + 


2  \  1/2 


a  =  b  11  +  ^  1  . 


(2)  The  above  approximation  for  the  radius  value  is  to  be  used  with 

2 

Naperian  trigonometry  (assumed  to  be  valid  up  to  e  correction). 

10.3  Problems  with  the  current  methodology 

a.  On -an  oblate  spheriod  the  radius  vector  magnitude  varies  from 

b  <  1  r  I  <  a. 


b.  Now  since  e,  the  eccentricity  can  be  defined  as  e 
then  rearranging 


uf/i  2, -1/2 

a  =  b(1  -  e  ) 
and  if  e  <<  1  then. 


as  b  1 


therefore: 
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b  < 


-»■ 

r 


c.  For  the  equal  area  compensation 


=  b 


TECHNICAL  NOTE 


This  value  of 


is  assumed  to  be  used  over  the  whole  spheriod. 


d.  A  problem  arises  in  the  level  of  accuracy  required  for  precision 
measurements.  The  so  called  equal  area  method  should  have  an  accuracy  limita- 

2 

tion  on  the  order  of  (or  no  less  than  an  order  of  magnitude  less  than)  e 
2 

(on  the  earth,  e  »  .0067).  For  quarter  turn  measurements  (the  distance  swept 
out  by  a  90°  central  angle,  this  to  be  redefined  later)  errors  from  -5  to  +11 
kilometers  can  be  obtained. 

10.4  Introduction  to  FROST 


a.  Start  by  rederiving  the  spherical  case  geodesics,  considering  each 
necessary  component  separately,  and  relating  each  component  via  a  parametric 


(1)  For  the  spherical  case: 


variable,  A. 

(1) 

p(X) 

=  a 

(X ) 

=  tan  ^ 

9(X) 

+1 

o 

CD 

II 

I|j  (X  ) 

=  tan  ^ 

■1  ^ 
3 

2 


\  +  1 


1/2 


(synthetic  law  of  sines) 


As  (1 )  =  a  tan 


■l( 

\  1  +  3^ 


/  X 

1/2  I  1 


(synthetic  law  of  cosines). 


(a)  It  can  be  shown  that  upon  elimination  of  X,  the  above  rela¬ 
tions  reduce  to  standard  spherical  trigonometry. 
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(bj  Spherical  case  quadrant  analysis  for  t 


QUADRANT  I 


QUADRANT  II 


niRECT 


RETROGRADE 


for  the  heading  or  direction  angle,  tan  ^  3 
the  following  terms: 


1 

O 

-3" 


1/2 


a(X) 


Quadrant  I 

Quadrant  II 

q(X ) 

■rr  -  q(X)  DIRECT 

IT  +  q  (X  ) 

2ti  -  q(X)  RETROGRADE 

(2)  For  the  slightly  oblate  spherical  case  assume  the  eccentricity 
value  e  is  small.  Then  using  calculus  of  variations  and  salving  for  the 

2  Lj. 

geodesic  curves,  keeping  terms  up  to  e  only  (eliminating  e  and  higher  terms) 
obtain  the  following  relationships: 


P  (X  )  =  b  1  +  -I 
\  2 


1  +  X‘ 


(X)  =  tan  ^  A. 


0(x)  =  ±  sec‘^  J 


^  ( 1  + 


6  . 


1  +  3 
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(A )  =  tan 


>  1  \ 
x^  -  j 


As  (A  )  =  b  <1  + 


1/2 


4(1  +  3^  ) 


tan 


1  +  3 


(x^ 


1/21 


^(1  +  3^  ) 


2  ^  iTY 


A  +  1 


(a)  (j)  in  spherical  coordinates  is  defined  as  the  geocentric 
colatitude  variable.  The  geodetic  latitude  variable  is  related  to  A  by 

L(A)  =  cot  ^  "  e^)  aJ  . 

(b)  In  the  above  equations,  in  the  limit  as  e  ->■  o,  the  spheri¬ 
cal  case  solutions  are  recuperated. 


b.  FROST  makes  a  first  order  attempt  at  solving  for  the  geodesics  on  a 
slightly  oblate  surface. 


(1)  It  assumes  the  surface  eccentricity  (e)  is  so  small  that  terms 

4 

of  the  order  of  e  and  higher  may  be  discarded*  This  seems  a  reasonable 
assumption  since: 


turn* 


(a)  e  corrections  are  from  -5  to  +11  kilometers  per  nuarter 


4 

(b)  e  corrections  are  on  the  order  of  tens  of  meters  per 
quarter  turn,  assumed  to  be  well  VA/ithin  the  noise  of  the  actual  surface  fluc¬ 
tuations  of  the  geoid;  therefore, 

4 

(c)  Higher  order  terms  (greater  than  e  ]  can  be  absolutely 
discarded  for  the  earth  for  any  reasonable  accuracy. 


(2J  The  above  assumptions  form  the  basis  for  the  limitations  on  the 
use  of  FROST. 


10.5  Comparison  study 

a.  In  comparison  with  the  spherical  relations  note  an  apparent  crunching 
up  of  the  curve  for  9  at  the  nodal  crossings. 


9  dim  A  00)  =  ±  I 


sphere 


9  dim  A  ->■  oo)  =  9  ±  L 

o  I 


+ 


oblate 


correction 


also,  i|)(node}  yields: 
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-1 


4i(lim  X  =  tan  3  sphere 

2  \ 


dim  X  ->•  ")  =  tan 


-1 


1  e 


NOTE:  ipdim^  >  "}  =  ^  in  both  cases. 


oblate. 


A  I  ■'f  U 

As  I  =  a  •  sphere 


,  |X=«>  .(l+(l"e 


1  + 


/  2  2  > 
A  |X=~„irb/^  e  3 


4(1  +  3^  ) 
2  \ 


oblate. 


1  -  ^  - 


4(1  +  3^ ), 


Ti  b  /  A  ^ 

2  4 


1  +  3' 


(keeping  up  to  e  terms  only]. 


(1)  At  first  glance,  this  may  look  like  nonsense,  as  3  gets  larger, 
the  distance  traveled  is  smaller,  but  look  at  9(X).  As  3  gets  larger  Ad  =  tt 
is  not  transverse  so  that  this  result  does  make  sense. 

(2)  Parametric  geodesics,  symmetric  about  an  arbitrary  0^. 

(a]  Sphere. 
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(b)  Slightly  oblate  spheroid: 


TECHNICAL  NOTE 


h 


A0  =  7T 


b.  Comparison  of  distance  formulas.  Oblate  spheroid,  perturbed  spheri¬ 
cal  routing  with  great  circle,  equal  area  routing.  Comparison  will  be  from 
the  high  point  in  the  trajectory  (or  maximum  latitude]  (  X  =  g)  to  the  equator 
(X  =  »)  so  that  is  evaluated  as  B  <  X  <  <». 


(1)  For  the  perturbed  solution. 


A0 


2 


1/2 


-  e 


1  +  e 


X=oo 
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(2)  For  the  spheroid  approximation: 


As  ,  IT  ...  IT 

cos  —  =  COS  cos  •=•  +  sin  9  sin  ^  cos  AO 

c3  ^  ^ 


As  =  a  cos  Msin  (j)  cos  A0^ 

2 


but  cos  A0  =  sin 


Tie 


4  (1  . 32)1/2  • 


It  is  convonisnt  to  call  0  in  these  cases  tan  9^,  the  minimun  Golatitiirie 

value  for  a  given  geodesic,  so  that: 

2 


As  =  a  cos 
a 


“1 


•  X  ’  -A 

sin  4)q  •  sin  I  —  sin 


A  U  ^  2  ,  \l/2 


TTa 
-  2 

2 


e  •  2  , 

—  sin  ())g 


e^  e“  2 
1  +  ^ - T~ 


irb 
Ol  =  2 


1  + 


(up  to  e  terms). 


Now,  with  e  =  .0067227, 


a  =  6371.227  km 
b  =  6356.912  km  . 
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(3)  For  the  quarter  turn  (defined  as  3  <  A  <  : 


<0o 

As 

n 

0 

10002.196 

10007.899 

ID 

10001.690 

10006.805 

20 

10000.233 

10003.964 

30 

9998.000 

9999.489 

45 

9993.805 

9991.079 

60 

9909.609 

9982.669 

75 

9986.530 

9976.513 

90 

9905.414 

9974.259 

All 

figures  in 

kilometers. 

As 

n 

=  As  at  <|)^ 
a  o 

=  35.566^^ 

Max 

difference  j 

As  “  As  1 
n  a 

positive  =  11.155  km 

As  -  As 
n  a 

negative  =  5.703  km 

The  approximation  is  good  to  two  significant  digits. 

c,  ThB  comparison  of  FROST  versus  spherical  trigonometry  is  shown  in 

figure  1 . 

10. f 

3  Geodesic 

measurement 

example 

a.  Measure  distance  between  two  theoretical  bases: 


(1) 

Base  A: 

4125000  N,  7321146  W 

(New  York}  geocentric  latitude. 

(2) 

Base  B: 

6020000  N,  2920000  E 

(Leningrad)  geocentric  latitude 

(3) 

Alongitude  =  10229146 

(a)  4126910  M  geodetic  latitude. 

(b)  6021670  M  geodetic  latitude. 
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b.  Comparison  of  distances  measured. 

2 

(1)  Spherical  trigonometry,  equal  area  e  compensation,  geocentric. 
As  =  6744.0840  km. 

2 

(2)  Spherical  trigonometry,  equal  area  e  compensation,  geodetic. 

As  =  6716.2297  km.  A  =  8.7km 

(3)  FROST  (0  =  .460749734,  9^  =  7^2318).  .  „ 

o  A  =  1 9 . 1 km 

As  =  6735.3776  km. 

2 

(a)  Spherical  trigonometry  w/e  compensation,  geocentric. 

Let, 

A  be  at  longitude  o® 

B  at  longitude  102° 9 146 
A  =  {.748955721,  0,  .662620048} 

B  =  (-.111749248,  .487352137,  .866025404} 

0  =  60^64952234 

As  =  b  (^1  +  )  (1.058533854) 

=  6744.0840  km. 

NOTE:  0  =  cos”^  (A  •  B) 

2 

(b)  Spherical  trigonometry  w/e  compensation,  geodetic. 

A  =  (.746732702,  0,  .665124253} 

B  =  (-.111184443,  .484888955,  .867479522} 

0  =  60;39902851 

As  =  b  (^1  +  (1.054161913) 

=  6716.2297  km. 
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NOTE:  0  =  CGs"^  (A  •  B) 

(c)  FROST  calculation  (9^^  =  7^2318,  g  =  .460749734] 

A,  =  1.130294,  =  1.277565 

B,  Xg  =  0.577350,  Xg  =  0.333333 


As 

a 


b  tan 


p. 277565  -  0.2122903 
V  1.2122903 


0.006722  \ 

411 .2122903]  J 


0.006722 
4  ( 1.2122903  ) 


(l. 277565  -  .2122903  ) 
2.277565 


-1  /. 333333  -  0.2122903  0.006722  \ 

AS|^  -  D  tan  1.2122903  /  V  ''  4(1.2122903)  J 


0.006722  .  (.33333  -  .2122903  ) 

4(1.2122903)^'^^  1.33333 

As^  =  b  {.754145026  -  .000691664}  =  b (. 753453361 ) 
As.  =  b  {.306480974  -  .0003982597}  =  b { .306082714] 


As.^  =  b(1. 059536075) ,  with  b  =  6356.912 


=  6735.3776  . 


Given, 

g  =  .460749734 

A  =  1.130294 
a 

=  0.577350 
b 


calculate  A9^ : 


-1  A  e  3 

A9  =  sec  “  “  ^=1 - 1...  .  -  i 

3  2  r  .  .  .,2  M/2 


(  1  +  3^ 


an 


-  3^ 

1  +  3^  / 
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A0|^  =  65.94351463  -  .06068861  S  65.882826 
A0|^  =  37.05647301  -  .024663507  =  37.031809 
A0^  =  102.9146(351  . 
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APPENDIX  II 


20.0  THE  TORUS 


The  torus  is  a  very  complicated,  multiple  connected  surface,  and,  as  a  result, 
one  should  not  expect  nice  closed  form  solutions. 

What  is  presented  here  is  just  an  outline  of  how  one  might  go  about  solving 
for  the  geodesics  on  such  a  surface.  The  solutions  presented  are  strictly 
formal. 

Define  the  surface  model: 


r  i2  2  ,  2 

[r  -  aj  +  z  =  b 


(a  >  b) . 


Then  let. 


1  + 


r  =  a  i 


(l  ♦x2)1/2 


_ bA _ 

(l  .  x") 


(l  *  )’ 


(1  *x2) 


,  r  =  1  ± 


2^/2  ' 


g  =  9  ,  g  = 


Since 
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/»2  ^ 
^  g 


■=  C 


e  -2  '2 

+  r  g 


(1 


1  +  X 
and 

2  2 

-4  '2  C  E 

^  g 


1/2 


g  =  ± 


(1 

ce  1 


2  2-2 
+  r  c  g 


2  -  2  \  1/2 


1  +  X  r  I  r  -  c 


(  |rl  >  Ic]} 


Continuing, 


As  =  a 


As  =  a 


2  2 
c  e 


1/2 


1/2 


dA 


r  -2  “1 

e 

r 

.  2 

-2  2 

+ 

>- 

1 

r  -  c 

dA  . 


Finally, 

■X  dz  '>  dr  d9 

D=*n^i  +-7^-1  +r-priQ 

dA  z  dA  r  dA  9 


5(x)  =  ± 


b 

y\ 

1  + 

bX 

"i  4- 

1  — 
z 

[(1  * 

1  ^ 
r 

;  ±  b(i .  x^r" 


1  +  X 


r  (  r  -  ) 


6(X)  =  ±  1  ±  Xi  ± 

z  r 


2  1  1/2 


ce 


(r^  ^ 


r  V  r  -  c 


•0  • 


This  can  be  simplified  somewhat  by  letting  X  =  tan  i|/,  then 
r  =  a  ±  b  cos  t 
z  =  ±  b  sin  ip 
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9(ii))  =  0 

As('|<)  = 

5(i|;)  =  ± 


+ 


dii^ 


-  /  -2  2 
rV  r  -  c 


) 


1  2 


1  ±  e  cos  ijj  ^ 


i  ±  tan  '|»i  ± 

z  r 


f' sec  i|> 

Lv  " 


-  c"  ) 
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APPENDIX  III 

30.0  A  STUDY  ON  SURFACE  CURVATURE 

Surface  curvature  may  be  evaluated  using  the  parametric  geodesics  and  the 
appropriate  Frenet-Serrat  equations  because  the  geodesics  for  a  given  surface 
contain  within  them  qualitative  as  well  as  quantitative  information  regarding 
the  surface  itself.  The  parametric  geodesics  form  a  natural  network  or  grid 
which  define  the  surface. 


This  evaluation  is  started  by  presenting  the  Frenet-Serrat  equation  for  cur¬ 
vature,  evaluating  the  Christoffel  symbol  of  the  second  kind,  showing  the 
parametric  evaluation  for  the  symbol,  then  proceeding  directly  to  the  spheri¬ 
cal  case. 


i  i  dx^ 

If  X  are  the  coordinate  variables,  X  =  the  unit  tangent,  then 


dX 

ds 


r  1  ^  dx^  i 

X  -; —  =  Kli  (where  K  is  the  curvature,  y'’"  is  a  unit  vector). 


ct,3 


ds 


Now, 

ml 

2,2 

r-l 

3  ,3 
1-2 
1,2 
2 
3,3 


=  -r  , 


.  2  . 
sin  ^  , 


[ 


r 


=  -sin  (()  cos  <j)  , 


r3 

=  cot  (j)  .  All  other  F  ^  HQ 

2,3  'a, 3 

For  the  sphere  (r  =  a  ]  and  , 

p  1 

=  -a  , 

2,2 
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3  ,3 
^2 
1,2 
r-2 

3  ,3 
I—  3 

1.3 
I—  3 

2.3 


=  -a 


1  +  A' 


1_ 

a 


1  +  X 


2  ' 


a 


± 

X 


30.1  Spherical  case 
Letting 


X  =  a 


2  ^  "1  ^ 

X  =  tan  X 

3  1  A 

X  =  9q  +  sec  (choosing  the  positive  side), 

3  “  1  A 

NOTE:  The  result  is  the  same  if  x  =  9^  -  sec  -g-  is  chcosen. 


Then, 

A^  =  0 

,2  1  dX 


1  +  X 
3  1 


2  ds 


_ ^ 

1/2  ds 


X  _  -  1 
\ 


and  with 
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6^  (  1  «■  ) 
.8^)x’  ■ 


Now 


Ky 


2 


8^  (  1  *  )  I 

aHl  •6^)x^  *  " 
aHl  *6")x' 


-  0  . 

dx3 

Evaluating  is  somewhat  more  difficult.  Remember  that  the  great  circle 
route  in  parametric  curvilinear  coordinates  is  a  multiple  defined  curve 
9  =  9q  ±  sec  ^  ^  ,  based  on  a  multivalved  function. 

Consider: 


The  first  derivative  is  ± 


d9 

dA 


,  depending  on  which  side  of  the  curve  (or. 


better  said,  which  branch)  is  chosen  for  usej  however,  the  second  derivative 

2-2-  is  both  equal  in  magnitude  and  direction  (all  concave  down)  on  either 
dA^  „ 

side,  so  that  a  straight  forward  brute  force  — -  evaluation  gives  contribu- 

dA^ 


tions  from  both  branches.  The  term  must  be  divided  by  two  for  use  in  the  cur¬ 
vature  equation. 

Continuing  with. 
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dX 

ds 


=  (-)- 

V2  /dX 


^  _ 3  ,  1  •*•  X^  \ _ 1 _  (l  f  X^)rx^  - 


-  (1\  e  (  2X^  -  2X  (  1  +  x^)\  (l  +  x^)(x^  -  3^) 

^3^)1  ^ 


i _ (1  ^  X^Xx^  -  3^)^/^ 


3^  (  1  +  3^  ) 


and 

Kw' 


.  (.1.J_^!Xll.-.  .6.',  ).*/.'  .  1  [xlj3]  .  ^  [,2,3- 


(1  .8^) 


i _ (i  >  x^Xx^  - 


(  1  +  3^  ) 


(x^  -  3^)^^^ 


.  0  . 1  ir_: _ ^ 

X  I  c  .  .2  1  1/2 


(1  ^  X 


1.  +  X' 


(1^3^)^/^  x= 


=  0  . 


/-1 


i  1  -  ^ 

So,  the  solution  for  Km  =  —  (  0  ),  constant  and  in  the  ~i  direction,  i'nis 

'^  \  0  /  •' 

is  known  to  be  generally  true,  not  only  for  the  great  circle  routes  but  for  the 
surface  of  a  sphere  [constant  curvature,  independent  of  where  the  observer  is 
(6q,  A)  or  in  what  direction  the  observer  is  looking  (3)  I.  This  is  an  irnpor- 

tant  result.  It  also  points  out  what  would  have  happened  had  not  1/2  of  the 

dx3  3 

term  been  taken  when  evaluating  the  curvature  (Ky  0  but  a  function  of 

A,  3  'which  would  not  make  sense). 


*=^or  completeness  lets  turn  our  attention  to  evaluating  the  torsion  of  great 
circle  routes. 

Fortunately  it  is  a  trivial  matter  to  show  that  t,  the  torsion  for  the  para¬ 
metric  spherical  geodesic  is  x  a  0. 

From 


dy 


ds 


i  3 

a  dx  /  i  IN 

y  —  =  -  (KX  *  TU  ; 


now 
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^  -  A® 
ds 


and 


What  is  needed  to  be  shown  is ; 

li“x^  +  KX^  =  0  =  TW^  . 
a3 

NOTE:  ail  always,  otherwise  p  =  0. 
For  i  =  1 , 


tp^  =  0  . 

For  i  =  2, 


TP^  =  0  . 

For  i  =  3, 
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=  0 

OED. 

For  the  minima  solutions  on  the  sphere  there  is  zero  torsion  and  constant  cur¬ 
vature,  independent  of  3  or  X.  It  is  probably  primarily  for  this  reason  that 
the  sphere  is  so  clean  to  work  on. 

The  end  of  this  section  will  concern  itself  with  the  evaluation  of  curvature 
on  a  slightly  oblate  planet.  Coordinates  to  be  used  are,  of  course,  spheri¬ 
cal. 


The  appropriate  expressions  are: 

1 


p  =  b  1  + 


2  ,  2 
e  X 


2  1  .X^ 


=  X 


“1  2 
=  tan  A  =  X 


Q  n  ””1  ^  S 

0=0^+  sec  tt  - 


^  tan-^^^ 


2  „2  \ 1/2 


'  Ml  ^  \l  .  3^ 


=  X 


Now, 


§-b(l  .(,  1 


2  -2 
e  X 


dX 
and. 


1  *  ‘  >■")(>■’-  - 


dp  be 

dX 


(1 


d({)  _  1 

1  +  X^ 


d9 

dX 


1/2  2 


X  _  -  1 

I  3^ 


(l  ^  1 


1/2 


(in  the  derivation  e 


4  + 


3‘ 


lerms  will  be  collected  and  ignored  when  justified). 


ds 


Cx^  -  3'  ) 


b(l  +(  1  -  e^)3^)  .  .  (  d  .6-^2 


1  .(1 
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,  1  i  2  (  -  6^  ) 

(1  *(1  -3^)6^)*' 


1  ♦  (i  ♦fi)x2 


ds 


a'  =^  = 


(i  ^  x^Xx^- 

b(l  -  (i  - 

_ e  ( 1  +  X 

b(l  ^(1  ‘  [l  .  (l  . 


1  dx 
X  “  ds 


' )  r(i  *  d£ 

.2  ~  2  ~  ds 


ri*r,  .4i)xn  ^ 


Now,  the  terms  need  to  be  formed,  remembering  to  take  ^  f ^ 


After  considerable  algebra,  and  keeping  terms  up  to  e  only: 


3^  (  1  +  X^  )^  (  1  -  X^  ->•  28^  ) 


b(l+(l-e^) 


3^)n.(i.|!)x^T 


3^  (  1  X^  ) 


ds  u2  ^  .  f  ^  2X2X3 

b  V1+Vl-e  ) B  J \ 


i).^r 


X^  (  1  ^  X^  133^  H.  X^  (  3^  -  3 

2bni  (l  -  e^)3^)x^"  ^  si 


„  f  ^  ,2X2/', 2  „2X1/2 

-3X1'*'^  /  W  ~3  / _ 


b"  ( 1  *  ( 


(1  *X^)^.^X^(: 


Finally  the  curvature  can  be  calculated: 


jx  1  r  1  o 

dX  X  cix  3  ,,  1 

-r -  +  X  X  =  Kp 

ds  o 

a  ,3 


Again,  after  considerable  algebra,  the  following  is  obtained: 


.1  1  e^  1  X 

Kp  =  -  •r  +  T—  o  - 


^  fi  x^  ^  (l  +  23^  -  x^) 


b  b  2 


+  x^  (i  +  e^)(i  +  x^) 
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= 


bH  1  +  )  X 


.  (3S^  2X^3^  -  2s}  ) 

(l  6^) 


Km  =  0  . 

Combining  all  terms,  the  curvature  can  be  written  as: 


K  u  ^  +  K  u  ^ 
0^0  1^1 


where  is  an  e^  perturbing  term. 


s 


-1 


2  2 
1  +  2B  -  X 


A 


b  \  2  (l  +  X^)  (l  +  6^)(l  +  X^  )i 

33^  +  2X^3^  -  2X^ 
b^(l  ^  3^)(l  ^  xOx 


iKi  =(  g.j(  Km")(’kmJ^ 


30.2  Conic  curvature 

p  =  aX 

^  •  ”1 
9  =  3 1  n  a 


9  ■  ®0  ■  ^  ^  ^  ^the  positive  branch  is  chosen  and  as  in  the  case  of 


the  sphere  take  one-half  the  value  of 


d^9 

ds^ 


dp 

dX 


=  a 


^  =  0 
dX 


^  _  J _ 

dX  "  aX  /  2, 2 


\  1/2  ,  /  2,  2  2  1/2 

a  X  _il  aXvotX  -  c  J 
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ds  _  a  a  X 
dT  =  ^^^272  '2 


r-c-F^  • 


Now, 


,  Z'  2,  2  2X1/2 

dj^_(^aX  -  c  ) 

aX 


ds 


^  =  0 

ds 


d9  ^  _ c 

ds 


2,2 

a  a  A 


dp  ^  1 


1  r.  Z'  2, 2  2  ^  -1/2  2,  Z^  2,  2  2  'l  1/2  1  (a^X^  - 

-^|aX  (ax  -  c/  aX-(aX  -  c  )  -'  aj-  ^ ^  ' 


^2  2, 2 

ds  a  X 


1 


[aV  -  (o.h^  -  c')]  .  t 


2,  2  2X1/2 

X  -  C  j 


,  2  /  2,  2  2X1/2 

ctX  V  a  X  -  C  J 

-2  Z'..2, 2  _2jl/2 


a  a  X 


,  2  /  2,  2  2  X  1/2 

O.X  \  a  X  -  e  J 

2 


.  ( aH*~  -  c‘ 


a  a  X 


.  2^  3 

a  U  /\ 


2 

d  (|)  , 


5  □ 


ds 


1  d^9 


2,  2  2 

c  a  X  -  c 


2,2  2.  3 

ds  a  a  X 


aaX 


f  2,2  2X1/2 

(  g  X  -  c  / 


2  3,4 
a  g  X 


Then  with: 


dX^ 

“1 

ds 

a,S 

Ku^ 

2 

c 

2, 
aa  X 

Kp^ 

=  0 

,g  dx  ,,  i 
*  d?  ■  ■ 


(+  0)  -  aa^X  ^ 


2  2 
c  c 


2  4,4  2,  3  2,  3 

a  a  X  aa  X  aa  X 


=  Q 
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Kw^  =  0  +  0  -  cx  ( 1  -  ^ 

2,2  2  \  111 


2  4,4  • 

a  a  X 


Kii^  = 


c(aV  -  ^  1 

aX 


2  3,4 
a  a  X 


aX 


2.2  2  \  l/2  (2.2  2  X  l/2 

c  ^  g  X  -c  )  ^  c  \  CL  \  -  c  ) 


2  3,4 
a  g  X 


=  0. 

So  that  for  the  cone, 


2  3,4 
a  g  X 


Km 


2  (  .  2X1/2/° 

1  _  c  g  V  1  -  g  ) 


2,2 

ag  X  , 


+  0 


2  4,4 
a  g  X 


dim  g  ^  1 ,  K  s  0)  . 


Evaluating 

.1  2  / ,  2X1/2 

1  I  _  J_  ,  c  g  1  -  g  ) 


Ky- 


g  X 


NOTE:  !k|~^  . 

P 

It  is  easy  to  show  that  the  curvature  of  a  helix  (the  parametric  cylindrical 
geodesic)  is  in  the  i^  direction  only. 

For  the  three  cases  (sphere,  cone  and  cylinder)  where  it  was  shown  that  the 
surface  constraint  was  written  of  the  form  u^  =  constant  (i  =  1 ,  2  or  3;  u,  a 

curvilinear  variable),  the  curvature  Ky^*  of  the  represented  surface  (assuming 
it  existed)  had  a  component  ONLY  in  the  e^  direction. 


Whether  or  not  this  is  a  general  property  is  beyond  the  scope  of  this  tech¬ 
nical  note. 
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APPENDIX  IV 

40.0  SURFACES  OF  REVOLUTION  AND  LONGITUDINAL  ASYMMETRY 

Calculation  of  geodesics  on  surfaces  of  revolution  and  first-order  pertur¬ 
bation  due  to  longitudinal  asynmnetry  (analysis  in  spherical  coordinates)  can 
be  accomplished  as  follows. 

Starting  with  the  spherical  surface  metric: 


2  /^d(i>  y  2 . 2  Ad9 

,  (  ^  I  .  p  sin  4>(  ^  ) 


assume  the  surface  can  be  written: 


=  p((j))  ,  [  =  p(X)]  ,  a  surface  of  revolution. 


Then  with  <()  =  tan  X  ,  9  =  g(X)  obtain: 


As  =  p(X) 


1 

\  P  9(!> 


1  +  X' 


1  +  X 


As  =  I  p  (X  )  ( 1  +  q(X  ) ) 


1  +  X‘ 


1  +  X 


Since  the  partial  of  the  integrand  with  respect  to  g  =  0,  then,  as  done 
before,  let 


c  =  p  (X  )  •  ( 1  +  q  (X  ) ) 


1  +  X' 


1  +  X‘ 


-1/2  ,2 


where  c  is  the  first  constant  of  integration. 
Solving  for  g  gives: 


X  •  p(X)  2 

x^ - 


1  +  q (X ) 


L  p 


9  (X)  =  9  ±  C  .r - 

0  X  •  p (X) 


1  +  q  (X  ) 

X^  -  —  (l  +  x^) 
2 
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0Q ,  the  second  constant  of  integration. 

For  the  sphere  p  =  constant, 
q  =  0, 


letting  ^  =  c 


•  [xt(l  . 


c"  -2 

Letting  - -  =  g 

1  -  c 


..  1 
“  X 


1/2 


9  =  0g  ±  sec  as  before. 


For  a  slightly  oblate  planet,  where  e^^  terms  are  significant,  but  not  e*^ 
+  terms,  proceed  as  follows: 


P  (A  ]  =  b 


1  + 


+  X 


(j)  (X )  =  tan  ^  X 

calculate  6 (X )  by  noting 

.  2 


q (X  ]  =  e 


1  +  X' 


SO 


1  -  e 


1  +  X' 


1/2 


X  —  -  1 

\ 


g  =  1  -  e 


2  X 


1/2 


1  +  X 


1  +  e 


1  +  X 


1  -  B 


1  +  X 


1/2 
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2(1  b(i+x^)^ 


n  n  .  )  -1  X  ^  -1  /  X^  -  3^ 

F  -  -  (TT^ryTTi  ItTT'] 


AUA _ 

FT^ 


4  -  -  1 


/V  LJ/V 

Y^( 


1  +  X‘ 


This  solves  explicitly  as: 


6(x)  =  6^  ±  isec-^  I  -  f- 


2(1  •  li  lx"  \  1  *  r 


x"  -  6"  V'" 


1  +  3‘ 


J  +  3^  ) 


2  X  1/2 


-l/x"-6"\‘/" 

^  TTF 


iHX j  and  As(X)  can  now  be  solved  for  by  the  usual  method: 


,,  .  -if,  /x"  .  1  r  e"  x"  B-  I  X' 

,HX)..an  •  '  ■“'TTF  “iTTx" 


As  =  b  1  -  e 


2  3 


2  \  1/2 


1  +  3  / 

X _ (il^l 


1+3^/  \  ^  4(1  +  3^) 


2  „2  \l/2 


1«x" 


tan-'fiFd 
♦  x"  ®  I  1  *  s" 


4 

/  .  2 

,2 

e 

''  -1 

-  3 

8 

1  +  X^  ' 

M  +  3^ 

1  *  S' 


1F)''l 


X"  - 


(1  *  i"  )  "  \  1  *  s"  /  (1  ♦  x")(l  *  s"  )  \  1  ♦  8 


X"  - 


,-i/x"-6"\F1|‘ 


(1.6")"  \i*e 
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Longitudinal  asyrrmetry  is  a  very  complex,  highly  nonlinear  problem,  solvable  pri¬ 
marily  via  successive  approximation.* 

Start  by  assuming 

P  =  a(1  +  eg)  where  e  <<1 

In  the  analysis  assume  the  primary  interest  is  in  working  on  a  small  surface 
portion,  or  sectoral  slice,  where  the  above  assumption  is  approximately  true^ 
also,  assume  working  in  the  low  to  mid  latitude  region  is  of  primary  interest, 

2 

thus  excluding  analysis  in  the  polar  region,  ie,  e  ~  0  and  e  <<  - -  for 

1  +  X^ 


all  X  in  the  domain  of  interest. 

On  the  earth  such  an  example  would  be  the  geodesic  from  St  Louis,  Missouri 
(elevation  ~  1000  feet)  to  Denver,  Colorado  (elevation  ~  6000  feet). 

Continuing : 


^  2  /  2.2/d0\' 

"  P  idlj  "  Hhx  j 


1/2 


dX 


p  dX  ~ 


and 


As  =  a  (1  +  eg) 


2  -2 
eg  + 


1 


\  1  +  X' 

(he  integrand  can  be  written: 


2  ,2  ^^11/2 

"  - ^  g 

1  +  X^ 


dX  . 


(1  +  eg) 


1 


1  +  X^ 


1  +  X' 


1/2 


dX  . 


which  can  be  further  approximated  as: 


(1  +  eg) 


1  +  X^ 


1  +  X 


1/2 


*Reference  Einstein,  Law  of  Gravitation,  pgs.  328-334;  Lass,  Vector  and  Tensor 
Analysis. 
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9  5 

In  this  example,  letting  the  integrand  be  6,  we  see  immediately  ^  0  , 

Therefore  we  must  use: 

d  36  -  n 

dX  3/  ■  9g 

and  collect  terms  in  e,  setting  up  for  a  successive  integral  approximation. 
Continuing, 


36 

3g" 


1 


1  +  ^  1  +  X^ 

,2 


1/2 


(1  +  eg)  • 


1  +  X^ 


(l  +  X^  )  ^  1  +  X‘ 


172  • 


r.  ^L.  95  f  36 

Rewriting  the  equation  as  —^r  =J  ^ 

(Since  as  oO  we  must  obtain  the  unperturbed  spherical  solution,  and  letting 

2  2  “ 

X  -2 


1  +  X' 


1  H-  X' 


g  =  e 


7/2 


1  +  X^- 

The  system  becomes : 


*  eg 


■J~[  dX  ) 


dX  +  c 


1/2 


1  +  X^ 


g  -  c 


1 


1/2 


e  •  <  [  1 


1/2 


1  +  X  /  1  +  X 

Start  by  setting  e  =  0,  solving  for  g  obtain: 
-1  X 


eg 


gg  =  %  * 


3 


Substituting  this  on  the  right  side  of  the  equation,  obtain: 


1/2 


1  +  X' 


1  +  X 


1/2 
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Q  +  -1  ^ 

c  0Q  ±  c  sec  + 


1/2 


1  +  /  (x^  +  1  )  /  X' 


1/2 


dX 


g  -  c  [  =  e(X) 


1  +  X 

where  e  (X ) 


e  ♦ 


1  +  3 


2  \  1/2 


1  * 


The  next  iteration  for  g  (=  g^l  can  now  be  set  up  and  performed  numerioally 
(note  e (X )  has  a  ±  term,  this  accounts  for  asymmetry  in  the  ig  direction). 


Now  obtain; 

X^  (  ( 1  -  c^  )  -  c^  )  =  0^  +  2c  £(X)  •  (  1 "  (e‘  ~  0) . 

continuing  the  process  of  successive  approximation : 


1/2 


xHx^(l  -c^)-c2) 


g  =  c  +  2c  s:(X) 


1  +  3 


2  \  1/2 


■1  .  ^^2  \  ,2  ,2 
1  +  A  \  A  '  p 


,  .|e(X)  .-i_ 

1  +  X 


1  +  3 

2  2 
X^  -  3 


2  \  1/2' 


1/2 


and  finally: 


g  =  0  ± 

*  0 


(1  +  £  •  r(X)) 


1/2 


^  3^ 


1/2 


dX 


-  1 


where  r(X  )  =  —  ^ 


1  +  3 


^  (i .  x2)nx2  -  3^i 


Q  .  -IX  ^  -1  /  X  -  3 

c  0^  ±  sec  -5-  +  tan  - - 

“  ®  ll*6' 


2  \  1/2 


NOTE:  3 


1  -  c 


Caution  must  be  used  when  determining  which  combination  of  ±’s  need  to  be 
employed,  and  on  which  side  of  0q.  Once  this  is  accomplished,  i|)(X)  and  As(X) 

can  be  calculated  in  a  straight  forward  manner. 
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PREFACE 

This  is  the  first  addendum  to  AFSPACECOM  Technical  Note,  Surface  Navigation 
and  Geodesy,  A  Paraaetric  Approach,  1  March  1988.  This  addendum  will  address 
two  issues  not  covered  in  the  aforementioned  publication. 

First,  a  practical  problem  involving  conic  geodesy  will  be  presented  and  the 
solution  discussed. 

Second,  the  solution  for  the  parametric  geodesic  of  a  paraboloid  (parabola  of 
revolution)  will  be  presented,  but  the  analysis  will  involve  the  use  of  para¬ 
bolic  confocal  coordinates. 

It  is  hoped  this  addendum  will  further  demonstrate  the  value  of  the  use  of 
parametric  analysis  in  solving  problems  of  surface  geodesy. 


THIS  TECHNICAL  NOTE  IS  THE  WORK  OF  DR.  K.  CARLTON-WIPPERN  AND  DOES  NOT 
NECESSARILY  REPRESENT  THE  OFFICIAL  VIEW  OF  THE  US  AIR  FORCE. 
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I.  AN  APPLICATION  OF  CONIC  GEODESY 

In  many  circumstances  the  surfaces  of  certain  physical  features  can  be 
approximately  modeled  (to  a  very  high  degree)  by  that  of  a  cone.  A  case  in 
point  would  be  Mount  Fuji  in  Japan,  which  will  be  assumed  to  be  in  the  shape 
of  a  perfect  cone.  Consider  the  following  problem: 

A  group  of  Shinto  monks  are  making  a  pilgramage  from  their  home  in  Kofu  to 
their  seaside  temple  at  Numazu.  On  the  route  their  master  has  selected,  they 
will  have  to  pass  across  the  base  of  Mount  Fuji,  from  the  northern  most  point 
to  the  southern  most  point.  The  leader  of  the  noble  group  of  pilgrim  monks 
has  been  charged  with  taking  the  route  about  or  over  Mount  Fuji  which  will 
traverse  the  minimum  distance.  It  can  immediately  be  seen  by  symmetry  that 
there  will  be  at  least  two  routes,  one  about  the  east  face  and  the  other  about 
the  west  face.  Not  withstanding  this  temporary  dilemma  (which  will  be  decided 
by  the  toss  of  a  coin)  the  following  elements  will  have  to  be  calculated/ 
determined: 

a.  The  minima  curve. 

b.  The  total  distance  traveled. 

c.  How  high  up  the  side  of  the  cone  will  the  monks  need  to  travel. 

d.  The  initial  angle  of  attack. 

Given,  Mount  Fuji  has  conic  parameters. 


So  that, 

p  =  (h^  +  T^y/^ 
0 

_  r 

“-h 

r  = - ^ 

(h^^Tr^yiA 
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Solution.  Reference  paragraph  IV,  Surface  Navigation  and 
Approach,  and  obtain  the  parametric  conic  geodesic: 


p  =  aX 

0  =  tan"^  a  (constant) 


A  diagram  of  the  problem. 


The  selected  geodesic  is  symmetric  about  the  HIGH  point 
chosen,  thus  it  can  be  concluded  that  occurs  at  9^ 


Returning  to  the  parametric  geodesics; 


.*•  - 1  I  TX 

from  0  =  6„  ±  sec  (  — 


r 


,  n  \  _  max 
sec  r  T5 - 


(from  p  =  aX) 


X 


max 


a 
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SIDE 

TOP 

in  the  route  to  be 
±2  ,  so  that 
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so  that 


this  fixes  the  parametric  curve  (6^  =  |  for  west  face  ascent,  ^  for  east  face 
ascent. ) 

Now,  p  occurs  at  I  0  -  0  I  =  0,  so  that 
imn  '  0 


From  this  obtain 
'■min  =  \  (?  )  • 

The  distance,  measured  from  ground  level,  traveled  up  the  side  of  the  mountain 
is 


The  total  distance  to  be  traveled  from  one  side  of  the  mountain  to  the  other 
can  be  evaluated  as  follows: 


X 


So  that,  for  both  ascent  and  descent  parts  of  the  geodesic: 


p2  ^  i/z 

4,  =  2a(X^  - 


x=(p  /a)=X 
^^0  niax 


X=(p  .  /a)=X  . 
^'^min  •'  min 


Now,  since 

X  =  ^ 
min  r  ’ 
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^max  ~  iT  ’ 


fp' 

r^p"' 

=  2a 

0 

-  _ 1  COS^  1 

^  Tir\ 

y 

y^a^  ' 

)_ 

As  =  2p  sin 

The  initial  angle  of  attack,  defined  as  the  angle  as  measured  off  from  the  top 
of  the  mountain,  measured  in  the  plane  whose  normal  is  parallel  to  the  conic 
surface  normal . 


Summarizing  the  solutions: 

Total  distance  traveled  =  2p^ 

Distance  of  ascent  up  the  mountain  side 

Initial  angle  of  attack  =  ^  (1  -  T). 
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II.  SOLUTION  OF  THE  PARAMETRIC  GEODESIC  FOR  A  PARABOLOID  (SURFACE  OF  A 
PARABOLA  OF  REVOLUTION) 

This  section  will  be  brief,  outlining  only  the  most  necessary  steps  required 
to  get  to  our  objective,  that  being,  defining  a  surface  constraint  and  develop 
the  five  equation  parametric  solution  for  the  geodesic  curves  (the  necessary 
surface  navigation  relations). 

For  simplicity  sake,  the  surface  defined  in  cylindrical  coordinates  is: 


where  a  is  an  arbitrary  parameter. 

Parabolic  confocal  coordinates  can  be  described  as  the  set  (^,  n,  0)  such  that 
0  is  the  same  longitudinal  variable  as  described  in  cylindrical  or  spherical 
coordinates , -then  using  p  from  spherical  coordinates  and  Cartesian  z, 

C  =  P  -  z 

n  =  p  +  z 

and  the  surface  constraint  equation  can  be  transformed  to: 
r^  =  a(2z  +  a) 

+  2az  +  a^ , 

z^  +  r^  =  p^ 

=  (z  +  a)^ 


p  =  z  +  a 


p  -  z  =  a  =  C 

hence,  the  surface  constraint  is  equivalent  to  setting  one  of  the  curvilinear 
variables  to  a  constant,  namely  C  =  a. 


Since  the  Jacobian  matrix  is  formed  to  show  curvil inearity  and  obtain  the 
scaling  factors,  it  is  convenient  to  work  with  spherical  coordinates  as  an 
intermedium: 


/p  sin 
P  sin 
\  p  cos 


0  cos  0 
0  sin  0 


p  sin  0 
p  sin  0 
z 


cos  0 
sin  0 
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P  -  z  =  ? 

P  +  z  =  n 

2z  =  n  -  ? 
so  that 

z  =  Kn  -  0 

Now 

p  =  Kn  +  C) 
so  that 


cos  0  =  -  =  ^  '  I  . 
p  n  +  4 

sin  *  =  (1.(1^  )')■'" 

=  (ti^  +  2ti?  +  +  2ti^  + 

_  2(Tig)V^ 

"  T1  +  e 

X  =  p  sin  0  cos  9 
=  Kn  +  C)2  cos  e 

=  cos  9 

A  similar  relation  exists  for  y,  so  that 
X  \  /  cos  9 

sin  9 

iCn  -  ?) 

The  Jacobian  matrix  can  be  formed  by  taking  the  appropriate  partial s: 


r 


cos  9 


cos  e  -(n£)‘/'  sin  9 


sin  9  cos  9 

0 


1  /n  y/^ 

'  U  j 
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It  is  now  a  simple  matter  to  show  that,  letting  I.  represent  the  vector  of  the 
ith  column  in  J, 


‘i  •  'j  =  [^j 


where 


6„  =  (0  i  #  J 
'J  11=  J) 


Thus,  the  coordinates  are  truly  curvilinear  and  the  scaling  factors  may  be 
obtained  as: 


The  distance  equation  becomes: 

(ds)2  =  ^  (n  +  C)(dti)*  +  Ti?(de)2 


Rewriting  parametrically,  with  ?  =  a 


3X  ^’’MdX 


dX  . 


There  are  a  variety  of  ways  which  can  be  chosen  to  solve  this  equation,  for 
example 


Ti  =  af 


s  =  a 


^(f  +  l)f'2  +  fg'^]  ^/\  dX 


and  since  the  integrand  is  not  dependent  on  g  explicitly,  then  letting  6  repre¬ 
sent  the  integrand. 
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set 


—  =  c 
ag' 

[^(f  +  l)f'2  + 

=  c^  [-^(f  +  l)f'^  +  fg'^] 
[f^  -  fc^]  g'2  =  ^{f  +  l)f'^ 


Let 


lJ-L  =  x^ 

f  -  c^ 

f  +  1  =  X^f  -  X^c 
X^c®  +  1 


f  = 


f'  = 


X^  -  1 
1 

(x^  -  ly 

2X 

(X^  -  1)^ 


[(: 


=  -2X(c^  +  1) 
(X^  -  1)^ 


l)2Xc2  -  2X(X==c2  +  1)J 
-  X^^c^'  -  l] 


Continue  with 


g  = 


0 


0 


=  9 


0 


^  c  f-2X(c^  1)  .  (X^  -  1) 

■  2j  (x^  -  1)2  X^c^  +  1 


XdX 


=  9 


0 


±  c(c2  +  1) 


XjdX _ 

l)(X2c^  +  1) 
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1 _ _  =  _1_  r  1  .  _i _ 1 

l(ax  +  1)  1  +  a  l_x  -  1  ax  +  Ij 


so  that  letting 


a  = 


X  = 


g  =  0  ±c  X^dXf  ;  . .  ]  -  (  -r-7 - 

J  W  -  i)  Vx^c^  +  1 

Solving  the  integral : 


X^  -  1  c^X®  +  1 

dX  =  -  f  dX 

X^  -  1  J  1  -  X^ 


X  =  tanh  u 
dX  =  sech^  udX 


r — -  dX  =  -  ftanh^  udu 

J  X^  -  1  I 

=  -  J  1  -  sech®  udu  (sech^  u 


=  -  (u  -  tanh  u) 
=  (X  -  tanh"^  X) 


c^X^  +  1 


cX  =  tan 


c*X®  =  tan*  0 


dX  =  -  sec*  0d9 
c 


obtain: 
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-[Ijtan  ^  0de 

-  9  -  1  do 

-  (tan  0  -  0)J  =  -  I  cX  -  tan'^  cX  |j 

g  =  9^  T  c  |^(X  -  tanh’^  X  -  X  -  1  tan'^  cX  J 

0  =  9  ±  l^c  tanh"^  X  -  tan"-^  cX] 

Recapping  the  first  two  coordinate  variables: 
$  =  a 

T,  = 


In  this  example,  one  constant  of  integration,  c,  appears  in  two  coordinate 
variables,  however,  the  curve  is  still  symmetric  about  an  arbitrary  9^,  as  one 

would  expect  from  a  surface  of  revolution. 

Continuing  for  ip, 

S  =  h  ^ ;  -f  h,  ^  s 

1  dX  1  3  dX  3 


tan  ^ 


D  •  e. 


D  •  e. 


which  resolves  as 
’X^'c^  +  1 


1  =  a 


-  1 


dX 


dX 


•2aX(c^  +  1) 

(X^  -  1)^ 

c(c==  +  1)X^ 


(X^  -  l)(X2c2  +  1) 
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xfcijLi 

-  1  / 


Finally,  for  distance  transversed: 


-  1 


4(X2c' 


X^c*  +  1 


i_  .  /x!ci_Li  ,  y  . 
+  1)  V  x^  -  1  ; 


4xMc^  +  ly 

(X2  -  1)* 


cHc*  +  D^X** 


x2  -  1  V  (X"  -  1)"  (x^c"  +  1)= 


xVc^  +  ly  \  V 


(x^  -  1)^ 


a(c2  +  1)  =  a(c2  +  1) 

J  (X*  -  1)2 


X^dX^ 

(1  -  X2)2 


AS 


sinh  2  tanh'^  X  -  2  tanh‘‘  X 


X  w/appropriate  change  of  limits, 
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PREFACE 

This  is  the  second  addendum  to  AFSPACECOM  Technical  Note,  Surface  Navigation 
and  Geodesy,  A  Parametric  Approach,  1  March  1988.  This  addendum  will  address 
two  classes  of  problems  not  previously  covered. 

First,  a  dicussion  of  the  surface  overtake  problem  is  presented.  This  is  a 
classical  naval  problem,  whereby  a  submarine  is  sighted  by  an  intercept  craft; 
the  submarine  submerges  and,  at  an  assumed  constant  velocity,  races  away  at 
some  random  angle  6  (measured  with  respect  to  the  initial  line  of  sight  of  the 
intercept  craft).  Assuming  the  velocity  of  the  intercept  craft  to  be  greater 
than  that  of  the  submarine,  what  is  desired  is  the  curve  the  intercept  craft 
must  take  to  eventually  overtake  and  catch  the  submarine  (this  is  a  well  known 
curve).  The  presented  problem,  solved  on  a  local  planar  surface,  also 
addresses  certain  quantities  of  operations  research  interest,  such  as  mean 
distance  to  intercept,  variance,  mean  time  to  intercept,  etc.  The  overtake 
curve  for  a  sperical  surface  is  then  solved  for  and  presented. 

Second,  the  problem  of  minimum  distance  navigation  on  composite  surfaces  will 
be  addressed;  specifically,  navigation  on  a  sphere  which  contains  a  spherical 
dent  on  its  surface  (an  approximation  of  the  surface  of  the  martian  moon, 
Phobus) . 


THIS  TECHNICAL  NOTE  IS  THE  WORK  OF  DR  K.  CARLTON-WIPPERN  AND  DOES  NOT 
NECESSARILY  REPRESENT  THE  OFFICIAL  VIEW  OF  THE  U.S.  AIR  FORCE. 
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I.  SURFACE  OVERTAKE  PROBLEM 

This  is  a  special  application  of  the  general  curves  of  pursuit  problem 
(for  example,  reference  Differential  Equations  with  Applications^).  In  this 
case,  for  a  given  surface,  an  interceptee  vehicle  is  sighted  by  the  intercep¬ 
tor.  It  is  assumed  the  interceptee  takes  a  minimum  distance  route  away  from 
where  he  was  sighted,  at  some  random  heading  angle. 

a.  Solution  on  the  plane. 


We  set  up  the  coordinate  system  as  follows,  at  time  t 


the  inter¬ 


ceptor  is  at  point  A,  the  interceptee  at  the  origin  (point  B),  at  a  distance  d 
from  the  interceptor.  The  interceptee  takes  a  linear  route  away  from  the 
origin  at  velocity  v^,  so  that  at  time  t  =  0  the  interceptee  is  at  a  distance 

d  from  the  origin  but  at  an  unknown  (heading)  location. 


The  equation  of  motion  for  the  interceptee  is  (in  polar  coordinates) 


0=9  (unknown) 

0 

“  *"0  '^0^  (where  r^  =  d). 

^Ritger  and  Rose,  Differential  Equations  with  Applications,  p.  94. 
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Assuming  the  interceptor  travels  at  constant  velocity  v^  (v^  >  v^^),  the  polar 
equation  of  motion  for  the  inteceptor  becomes 
s2  =  f2  +  (s®  =  v^) 

and  for  eventual  overtake  and  capture,  the  interceptor's  r(t)  component  must 
match  exactly  to  that  of  the  interceptee,  thus 

v|  =  v^  +  (r^  +  Vj,t)2  02 
1  . 


Letting, 


u  =  1  +  ^  t 


i=  1 


l/z 


obtain 
1  du 


a 


u 


=  d0 


and  solving,  this  becomes: 
a  I A0  I 


u  =  u  e 
0 


Substituting  back  for  r 
a|  A6  I 


r  =  re 
0 


or  parametrically,  the  solution  curve  for  the  interceptor  is 

r  =  r„  +  v„t 
0  0 


(choosing  0^^  =  0). 
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b.  Statistical  results  of  the  planar  solution. 

With  the  interceptee  heading  angle  (e^j)  defined  as  a  random  variable,  it  is  a 

reasonable  question  to  ask  what  are  the  mission  statistics  for  the  average 
mission  (eg,  average  distance  and  time  to  intercept,  and  variance  values). 
Borrowing  heavily  from  Parzen^  we  look  at 


for  the  solution  curve,  then  letting  6  be  a  random  variable 

de  =  1^ 

a  r 


where  r  <  r  <  r  e 
0  0 


2na 


Normalization  for  the 


1  =  c 


zira 


dr 

r 


assuming  an  eventual  overtake  no  later  than  at  0  =  2tt. 
distribution  on  r  becomes 


9 


c 


1 

Zva 


so  the  density  function  p(r)  is 
r 


1 

Zva 


1 

r 


r  <  r  <  r  e 
0  0 


2  7ra 


P(r)  =  ) 


0  otherwise 

Solving  for  the  average  r 

<r>  =J’rp(r)  dr  =  -  1) 

Solving  for  the  variance  of  r;  first 

<r">  =  J’r^p(r)  dr  =  ^  (e*""*  -  1) 

and  the  variance 

52  =  <r2>  .  <r>2 


becomes 

^Parzen,  Modern  Probability  Theory  and  Its  Applications,  pp.  308-316. 
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y^2  ^ 

t2  _  r  ft  1  ^  a.  2  ZTfCL  1  4.  1  >1 

5  ~  |_6  V.1  ■  nay  ■*■  ira  ®  ■  Tia  j  J 

With  the  average  radial  distance,  the  time  to  intercept  becomes: 

<f'>  =  )f  +  V  t  ,  t  =  —  [<r>  -  r„] 

0  0  6  ^  6  V*" 

0 

where  t^  is  the  average  intercept  time  (v^^,  r^^  and  <r>  known). 
With  tg,  the  average  total  distance  traveled  per  mission  becomes: 


(1  +  2ita)] 


If  (hence  1/a  »  1  ), 


[Note,  ordinarily  the  density  function  would  be  a  function  of  s,  the  total 

CX0 

distance  traveled,  not  r;  however,  for  r  =  r^e  r  and  s  are  linearly 
related: 


=  1  - 


^  ^  ®  1  .  ..»  /d8  ^  ® 

dr 

Now, 


a6 

r  =  r  e  , 


dr 

de 


ar  e 
0 


a0 


ar 


and 
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so 


ds  =  dr 


l/2 


With  s  and  r  linearly  related,  the  developed  statistics  do  not  change.] 


c.  Solution  on  the  sphere. 


In  this  case  it  is  assumed  the  interceptee  takes  a  great  circle  route  at  some 
unknown  heading  away  from  the  origin  (defined  here  as  in  the  last  example, 
the  point  where  the  interceptor  saw  the  interceptee).  Let's  assume  the  origin 
is  at  some  known  latitude  and  longitude  {L^,  L^} .  We  rotate  the  sphere  so 

that  the  origin  is  at  the  north  pole  of  the  primed  system  (we  do  this  using  a 
topocentric  rotation  matrix®). 


sin  L  cos 

a 

-sin  L„ 

0 

cos  L.  cos 

Q 


sin  L  sin  L 
a  0 

cos  L„ 

0 

cos  L,  sin  L„ 
a  0 


This  rotation  produces  a  figure  similar  to  the  one  shown  below 


E 


Using  one  more  rotation  about  the  primed  polar  axis 


-sin  a 

cos  a 
0 


0 

0 

1 


®Mueller,  Bate  and  White,  Fundamentals  of  Astrodynamics,  p.  79. 
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this  aligns  the  interceptor  A  along  the  new  (0  =  0)  axis  as  is  shown  below 


E' 


where  (for  a  sphere  of  radius  a)  the  distance  between  B  and  A  is  d  -  a  •  0^. 

If  B  is  at  the  origin  at  time  t  =  -  (^-^)  then  at  time  t  =  0,  B  will  be 

0 

at  some  random  heading,  a  distance  d  =  a  •  away  from  the  origin: 


E' 


In  this  coordinate  system,  a  great  circle  is  represented  by  6  -  0^^  [consider  a 
vehicle  at  the  real  north  pole,  any  great  circle  route  would  be  due  south, 
along  a  line  of  constant  longitude].  Thus  in  this  coordinate  system,  the 
interceptee ' s  equation  of  motion  is 

0  =  0 

V 

0  ^ 

(b  =  (j)  +  —  t 

^  ^0  a 

Again  assuming  the  interceptor  travels  at  constant  velocity  v^  (v^  >  v^),  the 
equation  of  motion  becomes 
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=  p®  +  +  p®  sin®0  0  (and  since  p  =  a) 

V®  =  a®0*  +  a*  sin*0  0* 


In  analogy  to  the  planar  case,  the  0  components  of  the  two  curves  must  match, 
hence: 


v^  =  +  a^  sin^  (0^^  ^  ^  )0‘ 


v^  -  V®  \ 

_1 - 1) 

a2  / 


=  Sin  (  0^  +  —  t  )  0 


d0  =  dt  •  CSC 
Letting 


(k  *  r  0 


.  V? 


1/2 


U  =  0^  +  -  t 
and 


l/2 


a  d0  =  CSC  u  du 

1  +  cos  u 


a0  =  c  -  In 


sin  u 


In 


sin  u 
1  +  cos  u 


a0  -  c 


sin  u  _  A  «6 
1  +  cos  u  - 


but 


sin  u  _ 

1  +  cos  u 

and  setting  0  = 


at  t  =  0  obtain: 


0=2  tan"^ 
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or,  parametrically 

0(t)  =  F  ^ 


0(t) 


- 

- 

In  tan 

0  +  —  t 

0  O 

-  In  tan 

0 

2 

T 

- 

with  (0,  9)  solved  for  in  the  rotated  (or  primed)  system,  the  solution  in  the 
original  (unrotated)  system  can  be  easily  generated,  this  by  creating  the 
primed  vector 


r 


=  < 


sin  0  cos  6 

sin  0  sin  0 

cos  0 


then  finding 

P  =  p 


Interceptor  latitude  and  longitude  data  is  available  directly  from 


cos  L 


cos 


P  -  S 


cos  L 


a 

sin  L 


sin  L, 


Mission  statistics  (averages  and  variances  for  a  randomly  headed  interceptee) 
can  be  calculated,  albeit  numerically,  as  before. 

d .  Direct  intercept  problem. 

This  is  the  standard  curve  of  pursuit  problem,  whereby  the  interceptor  sees 
the  interceptee,  and  pursues  the  interceptee  along  a  curve  such  that  the 
interceptor's  heading  is  constrained  to  be  simultaneous  with  the  line  of  sight 
to  the  interceptee. 

(1)  Solution  on  the  plane. 
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Construct  the  coordinate  system  as  follows: 


y 


t  =  0 


where  A  is  the  interceptee,  initially  at  the  origin,  and  B  is  the  interceptor 
initially  located  a  distance  a  from  the  origin  along  the  x  axis.  Starting  the 
clock  at  time  t  =0,  then  at  time  t  =  6t  we  assume  the  interceptee  has  trav¬ 
eled  some  distance  es^^  along  the  y  axis  at  velocity  v^  and,  similarly,  the 

interceptor  has  traveled  some  distance  6s,,  at  velocity  y. 

Refering  to  the  figure  below: 


y 


In  time  6t  the  interceptor  has  moved 


fi(s„) 

And  the  interceptee  has  moved  — - —  in  the  same  time. 

*0 
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Now,  distance  s^,  measured  from  point  B  on  the  curve  y  =  f(x),  is  equal  to 
y  -  X  tan  ^  =  y  -  x  and  equating  the  interceptor  to  interceptee  movement 

in  time  6t: 


Letting  «  =  and  assuming  a  <  1  (v^^  <  vj 


d*y 

X  — -  =  a 

dx® 


1  M  ^ 


dx 


i/z 


with  initial  conditions: 


a) 


0 
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y(x  =  a)  =  0 


it  is  relatively  easy  to  show  the  solution  curve  is: 


y  =  a 


1 

1  +  a 


A 


(2)  Approximate  solution  on  the  sphere.  As  a  first  order 
approximation,  we  will  make  the  following  assumptions: 


(a)  The  interceptee  is  located  on  the  equator  and  traveling 
due  east  along  the  great  circle  route  [there  is  no  loss  in  generality  to  this 
assumption,  providing  the  interceptee  is  confined  to  travel  on  a  great  circle 
route.  This  being  assumed  true,  then  a  judicious  choice  of  coordinates  (ie, 
rotating  the  sphere)  will  put  the  interceptee  on  the  equator  as  desired]. 


(b)  The  initial  position  of  the  interceptee  is  analogous  to 
that  of  the  planar  case,  in  that  the  initial  line  of  sight  of  the  inter¬ 
ceptor  is  perpendicular  to  the  line  of  travel  of  the  interceptee. 


*Ritger  and  Rose,  Differential  Equations  with  Applications,  p.  512. 
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(c)  The  initial  displacement  of  the  interceptor  from  the  inter- 
ceptee  is  small  so  that  certain  trigonometric  approximations  can  be  made  [eg, 
spherical  position  of  the  interceptor  {in  0  coordinates}  is  such  that 
sin  0  »  1,  this  being  true  when  the  interceptor  position  is  close  to  the 
equator;  eg,  0  =  90°  ±  20°].  Continuing  as  before,  at  some  time  6t,  the 
placement  of  interceptor  and  interceptee  will  be  as  follows: 


Now,  for  point  B 


1  +  sin® 


d0  •  a 


where  a  is  the  spherical  radius,  and'^easuring  s^  from  point  B 


•^  =  9  +  tan’^ 

cl 


cos  0  sin  0 


d0/ 


The  second  term  is  obtained  from  spherical  trigonometry 


cos  I  =  sin  0  cos  0 
sin  9  =  sin  i  sin  0 
cos  9  =  cos  i  sin  0 

+  sin  Z  cos  0  cos  0 
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tan  0  =  cos  0  tan  0 

and  from  Carton-Wippern*,  where: 


tan  0  =  sin  0 


d0  ‘ 


Equating  both  sides  at  time  6t: 


1  +  sin' 


d0 ;  _ 


/z  1  f  JQ 

d0  =  —  d  0  +  tan"^  (  sin  0  cos  0  ^ 


Now,  using  the  assumption  from  Id(2)(c)  and  by  assuming  that  the  term  cos  0 

remains  small  so  as  to  approximate  tan’^  u  s  u,  the  equation  of  the  pursuit 
curve  becomes : 


''o  .  /d0  VI  _  d  /d0 

77L  J  -d?L®  "“"Ha? 


Continuing  on  the  right  side  with: 

^  +  cos  0  -  sin  0  ^  s  cos  0 

d0'  d0  ^  cl02 


and  letting 


®  ~  <  1)  > 


a[l  +  =  cos  0  ^ 

a  sec  0  d0  = - - 7- 

[TVp^ 

a  In  I  sec  0  +  tan  0  |  =  In  |  p  +  (p^  +  1)^/^ 
Letting  p  =  0  at  0  =  0^^,  then 


“Carlton-Wippern,  Surface  Navigation  and  Geodesy,  A  Parametric  Approach, 
pp.  38-43. 
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p  +  (p2  +  1)1/2 


tan_m 
tan  (xj 


where 


X 


Solving  for  p,  the  differential  equation  for  0  is 

^  ‘  ^^'*3 

Thus,  the  integral  equation  for  9  as  a  function  of  4,  assuming  convergence, 
becomes: 


4 

'  0 

where,  for  the  interceptor,  9=9^,  4=4^  at  time  t  =  0. 

To  show  that  the  integral  will,  in  fact,  converge  it  is  necessary  to  change 
variables,  0  =  |  -  L,  and  use  L,  the  conventional  latitude  variable. 

Then 


(tan  X|j)"“  •  (tan  x)“  -  (tan  x^)**  •  (tan  x)'‘*_ 


d0 


and  the  integral  becomes 


where  as  the  interceptor  approaches  the  interceptee,  L  (or  U)  -*■  0.  In  this 
limit,  it  is  easy  to  show  that  the  above  integral  (in  fact,  the  entire  spheri¬ 
cal  analysis)  converges  to  the  planar  case,  on  the  local  topocentric  plane. 
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II  OVERVIEW  OF  MINIMUM  DISTANCE  NAVIGATION  ON  COMPOSITE  SURFACES 

In  this  section,  we  will  confine  our  attention  to  solving  for  the  curve  of 
minimum  distance  between  two  points  when  the  surface  on  which  the  solution  is 
sought  can  be  described  as  two  simple  surfaces  (eg,  a  sphere  cut  in  half,  so 
that  the  surface  is  a  hemisphere  bounded  on  top  by  a  plane).  We  let  S^  be  the 
first  simple  surface  with  geodesic  be  the  second  simple  surface  with 

geodesic  and  let  c  =  c(u,  X)  be  the  locus  of  points  which  make  up  the  bound¬ 
ary  between  the  two  surfaces.  If  the  geodesic  between  two  points  traverse 
the  boundary  c  at  say  (u^,  Xj,  then  rp^(u^,  Xj  =  ^^(u^,  Xj  by  continuity. 

Futhermore,  if  the  boundary  is  smooth  and  continuous. 


T 

1 


dip^ 

dT 


(u, 


K)  =  T, 


the  tangent  must  also  be  continuous.  While  these  conditions  are  necessary, 
they  are  not  sufficient.  The  problem  is  a  nonlinear  one,  requiring  judicious 
evaluation. 


Example.  Let  P^,  P^  and  P^  be  three  points  on  the  equator  of  a  sphere  equally 
spaced  (ie,  P^  at  0®  longitude,  P^  at  120°  longitude  and  P^  at  240° 
longitude).  The  geodesic  (ip^)  between  P^  and  P^  lies  along  the  equator,  from 
P  ,  towards  the  east.  Similarly,  the  geodesic  (i>^)  between  P^  and  lies 
along  the  equator,  from  P^  towards  the  east.  Patching  the  two  minima  curves 
tp^  and  tp^  together  at  P^  it  is  observed  that,  at  P^,  '#'i(P2)  =  '^’2(^2) 

T,(P3)  =  T^CP^). 


but  the  true  geodesic  from  P^  to  P^  is  not  ^2* 
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The  geodesic  from  to  starts  at  P^  and  goes  westward,  around  the  back  to 

P  .  The  solutions  for  geodesics  on  more  complicated  surfaces  (such  as  a  com¬ 
posite  surfaces)  can  become  very  involved  and  are  not  necessarily  unique.  In 
the  next  section  we  will  address  a  simple  yet  practical  problem,  solving  for 
the  curve  of  minimum  distance  on  a  sphere  with  a  spherical  dent  feature  on  its 
surface.  Such  a  representative  surface  does  exist  (approximately)  in  nature; 
an  example  would  be  the  martian  moon  Phobus. 

Ill  A  PHOBUS  PROTOTYPE  COMPOSITE  SURFACE  NAVIGATION  PROBLEM 

Consider  the  following  diagram: 


sphere  of  radius 


R 


Thus,  the  spherical  dent  in  the  outer  sphere  of  radius  R  is  a  removed  partial 
sphere  of  radius  a.  A  useful  coordinate  system  is  one  which  describes  the 
boundary  between  the  two  surfaces  as  a  constant. 
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The  outer  sphere  has  its  center  at  the  origin,  the  boundary  between  the  two 
simple  surfaces  at  The  removed  sphere  having  its  center  at  A  and 

corresponding  boundary  at  ^  Letting  the  outer  sphere  be  with  geode¬ 

sic  and  the  removed  (inner)  sphere  be  with  geodesic  five  situations 

can  arise  for  any  two  distinct  points  (P^  and  P^): 

a.  P  and  P.  may  both  reside  in  S  . 

0  1  1 

-►  -► 

b.  P  and  P  may  both  reside  in  S,,  or 

0  1  2 

c.  P^  and  P^  in  and  respectively  (or  vice  versa). 

d  &  e.  P  or  P  or  both  may  reside  on  the  boundary  between  S,  &  S  , 

0  1  ^  . 

* 

Irrespective  of  the  case,  the  boundary  between  and  not  being  smooth 
requires  us  to  consider  three  curves  or  paths  when  solving  for  the  curve 
of  minimum  distance  between  points  P^  and  P^,  these  paths  are  ip^,  tp^*  and  a 
third  path  along  the  boundary  between  the  two  surfaces  (at  (P  =  <p^,  call  it 

tp^).  Assuming  P^^,  P^  are  unit  vectors,  then  in  S^,  the  distance  between 
points  is  AS^  =  R  cos"^  (P^  •  PJ,  similarly  in  S^,  AS^  =  a  cos"  ^  (P^  •  PJ. 
Example  1. 

Let  P^  be  on  the  outer  sphere,  P^  be  in  the  dent.  Find  the  minimum  distance 

curve.  Let  P(0g,  X)  be  on  the  point  on  the  boundary  where  the  minimum 
distance  curve  passes  through.  The  inner  and  outer  spheres  are  assumed  to  be 

longitude  aligned,  that  is,  P(0^,  X)  h  P(0^,  X).  The  idea  is  to  solve  for  X 
that  minimizes  the  total  distance. 

AS(X)  =  R  cos-^  [_Pjj  •  P(0^,  X)J  +  a  cos"^  Q  P^  •  P(0^,  X)J 

and  set 

9(As)  _  0 
3X 


*which  are  great  circle  curves  on  the  respective  spheres,  and  S^. 
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In  general,  a  minima  curve  between  points  may  cross  the  boundary  twice,  at 

P((>,  X)  and  P(0,  u).  From  this,  there  are  several  test  cases  which  must  be 
simultaneously  checked,  solved  for  and  compared  to  find  the  true  minima. 

Example  2. 

Let  P^,  P^  be  on  the  outer  sphere  (S^).  Under  certain  conditions,  the 
minima  curve  may  go  into  S^,  ride  on  the  border  or  remain  totally  in  S^. 

Sample  geodesics 


For  each  case  a  test  distance  is  used  AS^  =  R  cos"^  (P^  •  P^)  provided  all 
along  0,  0  ^  0^ . 
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AS^  =  R  cos-‘[^  •  P(«g,  X)]  +  a  cos-^[p(0^,  X)  •  P(0^,  u)]  + 

R  cos-^  [p(0^,  u)  •  P^  ^ 
and 

3AS,  3  AS, 

_ i  = _ i  =  0 

3X  3u 

ASj  =  R  cos-^[^P^  •  P(0^,  X)]  +  R  sin  0^  •  [  X  -  u  ]  +  R  cos-^[p(0^,  u)  •  P^J 

3AS,  3AS, 

— ^  =  - i  =  0 

3X  3u 

And  so  forth.  Once  P(0,  X),  P(0,  u),  etc.  are  found  which  satisfy  the  minimum 

distance  criteria,  it  is  an  easy  matter  to  go  from  P^  to  P^  by  way  of  dividing 

up  the  curve  and  using  the  appropriate  parametric  geodesic  relation.®  For 
navigation  on  the  inside  sphere  (S^)  it  is  convenient  to  think  of  the  south 

pole  as  the  north  pole,  ie,  solve  the  portion  as  though  one  was 
navigating  at  or  near  the  south  pole  of  the  excavated  sphere,  but  report  the 

position  latitude  or  w  -  0  or  northern  hemisphere  latitude.  In  this  way  the 
concept  of  north  pole  remains  consistent,  but  the  curves  must  be  reported 
separately  with  respect  to  each  simple  surface  (inner  or  outer). 

In  summary,  while  not  an  easy  task,  it  is  certainly  possible  (if  not  triv¬ 
ial)  to  examine  on  a  case  by  case  basis,  potential  geodesics  for  the  minimum 
distance  between  points  on  this  simple  composite  surface,  and  then  solving  for 
the  absolute  minima.  Depending  on  the  positions  of  A  (the  excavated  or 
removed  sphere  center)  and  the  ratio  of  a  to  R,  it  may  become  possible  to  make 

certain  generalizations  (if  and  are  in  the  geodesic  curve  may  always 
be  0j,  and  never  include  or  0^  curves;  also,  if  the  removed  spherical  por¬ 
tion  is  less  than  a  hemisphere,  0^,  the  border  or  rim  geodesic  vanishes  as  a 
potential  solution)  thus,  each  problem  must  be  evaluated  on  its  own  merits. 

K.  Carlton-Wippern,  PhD. 

Operations  Research  Analyst 


®Carlton-Wippern,  op.  cit.,  pp.  38-43. 
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PREFACE 

This  is  the  third  addendum  to  AFSPACECOM  Technical  Note,  Surface  Navigation 
and  Geodesy  a  Parametric  Approach,  1  March  1988.  This  addendum  addresses  the 
derivation  of  the  commonly  used  navigation  curve  called  the  rhumb  line  or 
loxodrome. 

The  loxodrome  is  a  well  known  curve,  constructed  so  that  the  angle  between  any 
meridian  and  the  loxodrome  is  constant  throughout  the  curve. ^  Although  there 
are  many  ways  to  derive  this  curve,  the  derivation  presented  here  is  based  on 
analysis  provided  in  the  subject  technical  note^  and  on  well  known  results 
from  general  curvilinear  analysis.  One  further  property  of  the  loxodrome 
should  be  mentioned;  when  plotted  on  a  Mercator  projection  the  loxodrome  be¬ 
comes  a  straight  line. 


^The  curve  of  constant  compass  heading. 

^Carlton-Wippern ,  Surface  Navigation  and  Geodesy  a  Parametric  Approach,  pp  4, 
33-35. 
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I  INTRODUCTION  AND  GENERAL  DERIVATION 

Consider  u  ,  u  and  u  curvilinear  coordinates  with  Jacobian  scaling  factors 

12  3 

h^,  and  h^. 


Then 

(ds)^  =  (h^duj^  +  (h^duj*  +  , 


and  if  u  ,  u  and  u  can  be  written  as  a  function  of  X 

1*2  3 

^  du  ^  du  ,  du  ^ 

D  =  h  -n—  e  +  h  -rr-  e  +  h  -rr-  e  ,*  where  D  is  the  tangent  (or  heading) 

1  dX  1  2  dX  2  3  ax  3 

vector. 

Assume  the  surface  under  consideration  can  be  constructed  by  setting  one  of 
the  curvilinear  coordinate  variables  to  a  constant,  u^=  constant  (eg,  for  a 

sphere  p  =  a). 

The  heading  angle,  measured  off  one  of  the  remaining  principle  axes,  can  be 
determine  by 

du, 

~  - ^  “  h^  *  HIT  rule).  (1) 

h  _2  2  2 

2  dX 

By  assuming  tan  i)  a  constant,  equation  (1)  becomes  the  loxodrome  generating 
differential  equation  which  relates  how  u^  is  linked  with  u^.  This  equation 

will  be  used  throughout  this  addendum.  Equation  (1)  can,  in  addition  to 
generating  a  loxodrome  curve,  be  used  to  generate  a  general  Mercator  projec¬ 
tion  of  a  surface. 

II  APPLICATION  TO  THE  PLANAR  SURFACE 
With 


obtain  two  loxodrome  forms  on  the  plane,  both  by  setting  z  =  constant.  For 
the  first  case: 

*Carlton-Wippern,  op.cit,  p  4. 
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tan  tp  = 


-  dy 


dx 


which  solves  as 
y  =  +  tan  (x  -  xj 

For  the  second  case: 


+  ,  de 

tan  T/;  =  r  ^ 


which  solves  as 


coti/)(9-e  ) 

z:  Y*  e  '  O' 
0 


(a  straight  line). 


(an  exponential  spiral). 


The  differential  equation  for  the  second  case  may  appear  to  the  reader  to  be  a 
fami 1 iar  relation .  ^ 

The  first  case  solution  (a  straight  line)  is  coincidently  the  planar  surface 
geodesic.  A  similar  analogy  exists  for  the  cylinder  (where  helix  is  both 
cylindrical  geodesic  and  loxodrome). 

Ill  APPLICATION  TO  THE  CONIC  SURFACE 

With 


(the  surface  generated  by  setting  0  =  a,  a  constant),  then 


tan  ip  =  p  sin  a  ^  , 

which  solves  as 

cot0»sina«(0-0  ) 

P  Pg  e  '■O'' 

with  0  oriented  as  discussed  in  Surface  Navigation  and  Geodesy.^  Distance 
measurement  along  the  conic  loxodrome  is  accomplished  by 


=  /  [1  +  tan^  0^  dp 


^Thomas,  Calculus  and  Analytical  Geometry,  p.  372,  equation  7a. 
^Carlton-Wippern,  op.cit.  p  19. 
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by  direct  substitution;  then; 

p 

As  =  sec  ip  Ap 

Po 

Singularities  occur  when  ip  =  it  (n  ±  with  n  =  0,  ±1,  ±2  ...  at  which 
point  p  =  pj  and 

0 

As  =  p  sin  a  •  A0 

00 

Also  of  interest  is  when  ^  =  nir  with  n  =  0,  ±1,  ±2  ...  at  which  point  6  =  6^ 
and 

p 

As  =  Ap 

Po 

IV  APPLICATION  TO  THE  SPHERICAL  SURFACE 
a.  Spherical  loxodrome. 

With 

and  p  =  a,  the  loxodrome  is  generated  by 

.  d9  * 

-tan  ip  =  sin  0  ^  , 

from  which 

r  4^  =  -cot  0  .  (9  -  0J 
J  sin  0  I-  V  0^ 


*The  minus  sign  introduced  to  reflect  the  proper  orientation,  <p  measured  from 
true  north;  therefore,  d0  becomes  -(d0)  so  as  to  produce  a  left-handed  coor¬ 
dinate  system. 
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.  * 

In  tan  ^  =  -cot  ^  •  (9  -  9  ) 

^  .  o' 

*0 

and 

tan  I  =  tan  ^  ^-cot*.(S-9J  ^ 

For  distance  measurement  along  the  spherical  loxodrome,  by  direct  substitution: 
As  =  a  [1  +  tan^  d^ 

<l> 

As  =  a  sec  if  A0 

0 

^0 

Again  two  important  singularities  arise;  0  =  it  fn  ±  i  jwith  n  =  0,  ±1,  ±2 

.  .  .  here  0=0  ^  ' 

0 

and 

0 

As  =  a  sin  0  A0  , 

0 

and  0  =  nil  with  n  =  0,  ±1,  ±2  .  .  where  0  =  0^ 
and 

0 

As  =  a  •  A0 


b.  Position  and  velocity  vectors  measured  in  inertial  space. 
By 

0 

As  =  a  sec  0  A0 


4 


1  June  1989 


TECHNICAL  NOTE 


changes  in  0  can  be  determined,  then  with 

^  .  “I 

9  =  +  tan  0  •  In  tan  ^  -  In  tan  |  (0  decreasing) 

changes  in  0  can  be  determined.  Now 
sin  0  cos  [9  +  fi(t-T)] 
sin  0  sin  [9  +  fi(t-T)] 
cos  0 

where  the  term  52(t-T)  has  been  included  to  correct  for  the  sphere  rotating  at 
an  angular  velocity  n,  measured  from  some  inertial  body  centered  coordinate 
system. 

For  the  velocity  vector,  let 

cos  0  cos  9  0  -  sin  0  sin  9  9 

cos  0  sin  90+  sin  0  cos  9  9 

-  sin  0  0 

0  is  determined  by  As  =  -a  sec  0  A0  (for  0  decreasing)  -s  =  a  sec  00  and 

d  9 

letting  S  =  v  (a  constant),  then  -  ^  cos  0=0;  and  with  -  tan  0  =  sin  0  ^ 

then  using  the  chain  rule  -tan  00  =  sin  0  0,.  and  0  =  ^  fln^ 

Now,  by  direct  substitution: 

-cos  0  cos  0  cos  9  -  sin  0  sin  9 

sin  0  cos  9  -  cos  0  cos  0  sin  9 

sin  0  cos  0 

We  add  to  this  vector  the  effect  of  rotation  and  obtain: 


-cos  0  cos  0  cos  9  -  sin  0  sin  0 
sin  0  cos  9  -  cos  0  cos  0  sin  9 
sin  0  cos  0 
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where  0=0+  ft(t  -  t). 

NOTE:  R  •  R  =  0  on  the  static  spherical  surface,  where 


V  APPLICATION  TO  THE  OBLATE  SPHERICAL  SURFACE 


a.  Slightly  oblate  spheroidal  loxodrome. 

We  start  by  introducing  the  oblate  spheroidal  curvilinear  coordinate  system: 
X  \  /cosh  u  sin  v  cos 
=  ^  [  cosh  u  sin  v  sin 
sinh  u  cos  v 


This  is  a  curvilinear  system  which  generates  oblate  spheroids  of  eccentricity 
e  by  setting  u  to  a  constant,"  specifically  u  =  sech"^  e.  Since  u  is  set 
constant  we  will  need  only  find  ha  and  h^  for 


As  =  ^ 


dX 


and 


hf  =  cosh^  u  cos^  V  +  sinh®  u  sin®  v 
2 


h®  =  cosh®  u  sin®  v  , 


so  with  cosh  u  =  - 
e 


As  =  ^ 
e 


(1  -  e®  si 


(e)’ 


+  sin®  V 


^  ^® 
dX 


i/z 


dX  . 


Now,  recognizing  ^  =  ae  (a,  the  semimajor  axis  of  the  ellipsoid)  and  con¬ 
tinuing  as  before: 


■tan  tl>  = 


sin  V 


d0 


(1  -  e®  sin®  V  *  dv 


"Which  is  equivalent  to  setting  curvilinear  (u^)  constant. 
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(I  -  e^^sin^  s,)±_  ^  ,  (3  .  B j 


Sin  V 


tanh 


- 1 


cos  V 


(1  -  sin^  v)^/^ 


-  e  tanh 


- 1 


e  cos  V 


(1 


s 


os  V  \ 

in^  ) 


=  -cot  i>  (9  -  0J 
and  with  «  1, 


— T-^ - S-  sin  V  dv  5  -cot  </»»(6-0) 

sin  V  2  '  O' 


V  G 

In  tan  2  Y~  ^ 


=  -  cot  )/)  •  (0  -  9  ) 


(for  V  decreasing) 


Thus,  this  is  the  first  order  correction  for  v  [which  by  the  way  is  not  <()  in 
spherical  coordinates,  but  is  related  by 

(1  -  cot  V  =  cot  0]. 

Distance  measured  along  the  loxodrome  becomes: 

As  =  a  J"  [1  +  tan^  •  (1  -  e®  sin^  v)^/^  dv 

=  a  sec  ij)  J'  (1  -  sin^  v)^/*  dv,  and  with  the  integrand  being 

g2 

approximately  1-2“ 


As  s  a  sec  ip 


Av  - 


Av 


2  sin 


2v 


)] 


Solution  singularities  and  resolutions  are  similar  to  those  on  the  sphere. 

b.  Position  and  velocity  vectors  measured  in  inertial  space. 

With  As  determining  changes  in  v,  which  then  determines  the  changes  in  0,  each 
successive  v,  0  is  found.  The  inertial  position  is  determined  by: 


a  sin  V  cos  [9  +  fi(t  - 
a  sin  V  sin  [9  +  52(t  - 
b  cos  V 


t)] 

0] 


where  a,  b  are  the  semimajor  and  semiminor  axis,  respectively. 
For  the  velocity  vector,  let 
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/  X  \ 

/ 

d 

dt 

y 

\zl 

'  \ 

a  cos  V  cos  0  V  -  a  sin  v  sin  0  9  \ 
a  cos  V  sin  0  v  +  a  sin  v  cos  0  0 
-  b  sin  V  V 

As  in  the  spherical  case  v  is  determined  by 
ds  =  -a  sec  (1  -  e®  sin^  v)^/^  dv 

s  =  -a  sec  (1  -  e*  sin*  v)^/*  v  , 
and  letting  as  before  i  =  v,  then: 

V  =  --  cos  ^  (1  -  e*  sin*  v)"^/* 

cl 

0  =  V  ,  sinji 
a  sin  V  ’ 

and  by  direct  substitution 

-cos  V  cos  0  cos  (1  -  e*  sin*  v)'^/*  -  sin  0  sin  ili 


dt 


I 


=  V 


-cos  V  sin  0  cos  </)  (1  -  e*  sin*  v)"^/*  +  cos  0  sin  ip 
(1  -  e*)‘/*  •  sin  V  cos  (1  -  e*  sin*  v)"^/* 


and  finally,  taking  into  account  rotation  about  the  polar  (or  minor)  axis; 


/ 


y  =  V 


-cos  V  cos  0  cos  (1  -  e*  sin*  v)’^/*  -  sin  0  sin  ip 
-cos  V  sin  0  cos  (1  -  e*  sin*  v)"^/*  +  cos  9  sin  4> 

\  (1  -  e*)^/*  sin  V  cos  ip  (1  -  sin*  v)"‘/*  J 


I  -sin  V  sin  9 
+  afi  I  sin  V  cos  9 

0  / 


where  0=0+  n(t  -  x) 
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VI  SURFACES  WHERE  THE  LOXODROME  AND  GEODESIC  CURVES  ARE  COINCIDENT 
Starting  with 

(ds)2  =  +  hiidu^y  + 

and  then  setting  u  =  constant,  thereby  generating  a  surface  spanned  by 

(u,.  U3); 


and  assuming 


ah^  /  aCh^  ;  h^) 

and  -  =  0  written  — - 7—  =  0 

3g'  \  9(g  :  g  ) 


the  equation  which  generates  the  geodesic  becomes: 


But  by  definition,  the  loxodrome  is  generated  by 
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tan  il>  = 


=  constant. 


Solving  both  equations  simultaneously  by  direct  substitution  yields: 


h 


tan  j) 
3  sec  rp 


c 


9 


h  =  c  •  CSC  ii>  , 

and  with  c,  tp  constant,  this  requires  that  h^  =  constant. 

The  geodesic  loxodrome  curve  relation  between  u^  and  u^  then  becomes: 

u,  =  u“  +  c"^  sin  tan  ^  I  h,  du,  , 

where  u“  and  c"^  are  the  two  constants  of  integration  from  the  (geodesic) 
second  order  differential  equation. 

There  are  surfaces  where  the  geodesic  loxodrome  exists  (eg,  on  the  plane, 
the  geodesic  loxodrome  is  a  straight  line;  on  the  cylinder,  the  geodesic 
loxodrome  is  a  helix).  This  mathematical  phenomena  can  occur  when  by  setting 
u^  =  constant,  to  generate  the  surface,  in  turn  implies  that: 

h  =  constant 

3 

and 


^  ^3) 

9(g  :  g') 


=  0 


This  is  not  to  say  that  the  proposed  process  of  generating  the  geodesic 
loxodrome  is  necessarily  unique  [3(h2  :  h^)/3(g  :  g')  =  0  may,  in  fact,  be 

over  restrictive],  but  that  further  analysis  into  the  phenomena  of  the  geode¬ 
sic  loxodrome  is  beyond  the  scope  of  this  addendum. 

On  a  variety  of  surfaces  (eg,  a  sphere)  a  parametric  geodesic  for  a  specific 
case  may  also  be  a  loxodrome.  For  example,  consider  the  geodesic  between 
two  points  on  a  sphere;  one  point  on  the  equator,  the  other  at  the  north  pole. 
The  great  circle  connecting  the  two  points  is  one  of  constant  heading  or 
direction  (due  north).  In  this  case  the  constant  c  =  0,  forcing  tp  =  0  (or 
180°),  and  the  geodesic  loxodrome  equation  reduces  to: 


where  0  is  the  longitudinal  variable. 
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VII  SUMMARY  AND  THE  TORUS  REVISITED 

The  loxodrome,  while  not  as  efficient  a  curve  as  a  geodesic  for  navigation 
purposes,  is  a  useful  navigation  tool  and  is  relatively  easy  to  construct  on  a 
wide  variety  of  surfaces  (provided  there  exists  a  set  of  curvilinear  coordi¬ 
nates  such  that  a  surface  under  investigation  can  be  generated  by  setting  one 
of  the  coordinates  equal  to  some  constant  value). 

As  a  final  comparison  between  the  tractability  of  solving  for  a  geodesic  or  a 
loxodrome  curve,  we  revisit  the  torus.  Choosing  curvilinear  coordinates  such 
that 


Figure  1.  Cross  section  of  a  torus 

In  this  coordinate  system  h^  =  a,  h^  =  av,  h^  =  a(l  +  v  sin  0),  and  setting  v 
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then  solving  for  9  =  9(0), 


d0 


®  j  (1  +  e  sin  0)[(1  +  €  sin  0)* 

However,  solving  for  the  loxodrotne: 


_  1  +  6  sin  0  d9 

•tan  0  = - -  TT 

e  d0 


-  c 


1+6  sin 


=  -cot  0  d9 


2e 


(1  - 


•  tan 


- 1 


tan  J  j  +  € 

L(1  -  J 


=  -cot  0  (9  -  9J 


0, 


NOTE:  The  torroidal  loxodrome  is  not  only  solvable  for  9  [  ie, 

2e 


9  =  9  -  tan  0 

0 


*  (1  - 


•  tan" 


tan  [  I )  ■*■  e 


L(1  - 


0, 


but  also  invertable  so  as  to  make  0  =  0(9);  given  0^,  9^^  ie,  calling 


2e 


(1  -  6^)V" 


tan" 


^an  (^)+  e 
_(i  - 


=  2  tan“ M  (1  -  •  tan 


g  -  e^)V^ 


2e 


(0„  -  cot  0  (9  -  0J) 


-  6 


and  the  measurement  of  distance  along  the  torroidal  loxodrome  becomes: 

l/2 


As  =  a 


*  (1  *  .  sin  ♦jM  s'' 


d0 


As  =  b 


1  ,  /1  +  €  sin  0  de 

^  - ; - Hi;  J 


l/2 


d0 


^remembering  e 


As  =  b  sec  0  A0 
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